ON THE SELF-SIMILAR ASYMPTOTICS FOR GENERALIZED 
NON-LINEAR KINETIC MAXWELL MODELS 
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A.V. BOBYLEV(*), C. CERCIGNANI(** ) , I.M. GAMBA(***) 

Abstract. Maxwell models for nonlinear kinetic equations have many applications 
in physics, dynamics of granular gases, economy, etc. In the present manuscript 
we consider such models from a very general point of view, including those with 
arbitrary polynomial non-linearities and in any dimension space. It is shown that 
the whole class of generalized Maxwell models satisfies properties which one of them 
can be interpreted as an operator generalization of usual Lipschitz conditions. This 
property allows to describe in detail a behavior of solutions to the corresponding 
initial value problem. In particular, we prove in the most general case an existence of 
self similar solutions and study the convergence, in the sense of probability measures, 
of dynamically scaled solutions to the Cauchy problem to those self-similar solutions, 
as time goes to infinity. The properties of these self-similar solutions, leading to 
non classical equilibrium stable states, are studied in detail. We apply the results to 
three different specific problems related to the Boltzmann equation (with elastic and 
inelastic interactions) and show that all physically relevant properties of solutions 
follow directly from the general theory developed in this paper. 
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1. Introduction 

The classical (elastic) Boltzmann equation with the Maxwell-type interactions is 
well-studied in literature (see [4,14] and references therein). Roughly speaking, this 
is a mathematical model of a rarefied gas with binary collisions such that the collision 
frequency is independent of the velocities of colliding particles. 

Maxwell models of granular gases were introduced relatively recently in [6] (see 
also [2] for the one dimensional case). Soon after, these models became very popular 
among people studying granular gases (see, for example, the book [13] and references 
therein). There are two obvious reasons for this fact. First, the inelastic Maxwell- 
Boltzmann equation can be essentially simplified by the Fourier transform similarly 
to the elastic one [5, 6] and second, solutions to the spatially homogeneous inelastic 
Maxwell-Boltzmann equation have a non-trivial self-similar asymptotics, and, in ad- 
dition, the corresponding self-similar solution has a power-like tail for large velocities. 
The latter property was conjectured in [16] and later proved in [8,9] (see also [3]). It 
is remarkable that such an asymptotics is absent in the elastic case (roughly speaking, 
the elastic Boltzmann equation has too many conservation laws). On the other hand, 
the self-similar asymptotics was proved in the elastic case for initial data with infinite 
energy [7] by using other mathematical tools compared to [8]. The recently published 
exact self-similar solutions [10] for elastic Maxwell mixtures (also with power-like 
tails) definitely suggest the self-similar asymptotics for such elastic systems. Finally 
we mention recent publications [1,15,22] , where one dimensional Maxwell-type mod- 
els were introduced for applications to economy models and again the self-similar 
asymptotics and power-like tail were found. 

Thus all the above discussed models describe qualitatively different processes in 
physics or economy, however their solutions have a lot in common from the mathe- 
matical point of view. It is also clear that some further generalizations are possible: 
one can, for example, include in the model multiple (not just binary) interactions 
still assuming the constant (Maxwell-type) rate of interactions. Will the multi-linear 
models have similar properties? The answer to this question is affirmative, as we shall 
see below. It becomes clear that there must be some general mathematical properties 
of Maxwell models, which, in turn, can explain properties of any particular model. 
Essentially, there must be just one main theorem, from which one can deduce all the 
above discussed facts and their possible generalizations. The goal of this paper is to 
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consider Maxwell models from a very general point of view and to establish their key 
properties that lead to the self-similar asymptotics. 

The paper is organized as follows. We introduce in Section 2 three specific Maxwell 
models of the Boltzmann equation: (A) classical (elastic) Boltzmann equation; (B) 
the model (A) in the presence of a thermostat; (C) inelastic Boltzmann equation. 
Then, in Section 3, we perform the Fourier transform and introduce an equation that 
includes all the three models as particular cases. A further generalization is done 
in Section 4, where the concept of generalized multi-linear Maxwell model (in the 
Fourier space) is introduced. Such models and their generalizations are studied in 
detail in Sections 5-10. The concept of an L-Lipschitz nonlinear operator, one of 
the most important for our approach, is explained in Section 4 (Definition 14. lj) . It 
is proved (Theorem 14.2)1 that all multi-linear Maxwell models satisfy the L-Lipschitz 
condition. This property of the models constitutes a basis for the general theory. 

The existence and uniqueness of solutions to the initial value problem is proved 
in Section 5 (Theorem 5.2). Then we study in Section 6 the large time asymptotics 
under very general conditions that are fulfilled, in particular, for all our models. It is 
shown that the L-Lipschitz condition leads to self-similar asymptotics provided the 
corresponding self-similar solution does exist. The existence and uniqueness of self- 
similar solutions is proved in Section 7 (Theorem 17. 1)1 . This result can be considered, 
to some extent, as the main theorem for general Maxwell-type models. Then, in 
Section 8, we go back to the multi-linear models of Section 4 and study more specific 
properties of their self-similar solution. We explain in Section 9 how to use our theory 
for applications to any specific model: it is shown that the results can be expressed 
in terms of just one function fi{p),p > 0, that depends on the spectral properties of 
the specific model. 

General properties (positivity, power-like tails, weak convergence of probability 
measures, etc.) of the self-similar solutions are studied in Section 10. This study also 
includes the case of one dimensional models, where the Laplace (instead of Fourier) 
transform is used. 

Finally, in Section 11, we establish in the unified statement (Theorem 11.1) the 
main properties of Maxwell models (A),(B) and (C) of the Boltzmann equation. 
This result is, in particular, an essential improvement of earlier results of [7] for the 
model (A) and quite new for the model (B). Applications to one dimensional models 
are also briefly discussed at the end of Section 11. 
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2. Maxwell models of the Boltzmann equation 

We consider a spatially homogeneous rarefied d- dimensional gas (d = 2, 3, . . .) of 
particles having a unit mass. Let f(v, t), where v G M. d and t G R+ denote respectively 
the velocity and time variables, be a one-particle distribution function with the usual 
normalization 

/ dvf(v,t) = l. (2.1) 

JE. d 

Then f(v,t) has an obvious meaning of a time- dependent probability density in M. d . 
We assume that the collision frequency is independent of the velocities of colliding 
particles (Maxwell-type interactions) and consider three different physical models 
(A), (B) and (C) described below. 

(A) Classical Maxwell gas (elastic collisions). In this case f(v,t) satisfies the 
usual Boltzmann equation 

/* = Q(fJ) = I dwdug f^f] \f(v')f(w') - f(v)f(w)} , (2.2) 



l xS d-l V \ U \ 

where the exchange of the velocities after a collision are given by 

v ' = ~(v + w + \u\cu) , and w' = -(v + w — \u\cu) 

where u — v — w is the relative velocity and Q G S^ 1 . For the sake of brevity we 
shall consider below the model non-negative collision kernels g(s) such that g(s) is 
integrable on [—1, 1]. The argument t of f(v , t) and similar functions is often omitted 
below (as in Eq. (JOJ))- 

(B) Elastic model with a thermostat This case corresponds to model (A) in 
the presence of a thermostat that consists of Maxwell particles with mass m > 
having the Maxwellian distribution 

'27iT\- d/2 ( mid 2 ' 
with a constant temperature T > 0. Then the evolution equation for f(v,t) becomes 



ft=Q(f,f) + » jdw<kjg{^\ [f(v')M(w') - /(»)%)] , (2.4) 

where 9 > is a coupling constant, and the precollision velocities is now 

v + m{w + \u\uj) v + mw — \u\u 

v = , and w = , 

1 + m 1 + m 

with u = v — w the relative velocity and uj G S^ 1 . 
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Equation (|2.4jl was derived in [10] as a certain limiting case of a binary mixture of 
weakly interacting Maxwell gases. 

(C) Maxwell model for inelastic particles. We consider this model in the 
form given in [8]. Then the inelastic Boltzmann equation in the weak form reads 

§-(f,^) = [ dvdwdwf{v)f{w)^j^[^{v')-^{v)] , (2.5) 

Ol JR d xR d xS d - 1 \ u \ 

where ip(v) is a bounded and continuous test function, 

f 1 + e 

(/, ip) = dv f(v, t)ip(v), u = v — w, lu G S d ~ l , v 1 = v — {u ■ lu)lu , (2.6) 

JR d 2 

the constant parameter < e < 1 denotes the restitution coefficient. Note that the 
model (C with e = 1 is equivalent to the model (A) with some kernel g(s). 

All three models can be simplified (in the mathematical sense) by taking the Fourier 
transform. We denote 

/(*,*) = ^[/] = (/,e-* w ) ) keR d , (2.7) 

and obtain (by using the same trick as in [5] for the model (A)) for all three models 
the following equations: 

(A) f t = Q(fJ) = Jdug(^) [/>+)/>-) -/(fc)/(0)] , 
where k± = \{k ± |jfc|o;), to G S d ~\ /(0) = 1. 

(B) f t = Q(f, f) +ejdug(^) \f(k+)M(k.) - f(k)M(0)] 



(21 



(2.9) 



Tlfch 



where M(k) = e~^t , k + = k. = k - k+, uo G S d ~\ /(0) = 1. 



(C) ft = / dJ-^[f(k + )f(k-)-f(k)f(0)} , 

(2.10) 

where k + = ^r(k ■ uj)cu, = k — k + , u G S^ 1 , /(0) = 1. An equivalent formulation 
is given for k=^-(k — \k\u), and k + = k — k_ . 

The case (B) can be obviously simplified by the substitution 

T\k\ 



/(M) = /(M)exp 
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Then we obtain, omitting tildes in the final equation: 



(B') 



ft = Q(fJ) + 9 dcog 



S d- 



k ■ uj 



f 



k + m\k\uj 



m 



/(*) 



(2.12) 

i.e., the case (B) with T = 0. Therefore, we shall consider below just the case (B'), 
assuming nevertheless that f(k, t) in Eq. (|2.12|) is the Fourier transform (|2.7|) of a 
probability density f(v,i). 



3. Isotropic Maxwell model in the Fourier representation 



We shall see that the three models (A), (B) and (C) admit a class of isotropic 
solutions with distribution functions / = f(\v\,t). Then, according to (|2.7|) we look 
for solutions / = f(\k\,t) to the corresponding isotropic Fourier transformed problem, 
given by 



x 



\k\ 



<p{x,t)=f(\k\,t)=F\f(\v\,t)] 



(3.1) 



where ip(x,t) solves the following initial value problem 



ft 



ds G(s) {<p[a(s)x](p[b(s)x] — f(x)} + 



+ 



(3.2) 



dsH(s) [tp[c(s)x] — (f(x)} , 

<Pt=o = fo(x) , f(0,t) = 1 , 
where a(s), b(s), c(s) are non-negative continuous functions on [0,1], whereas G(s) 
and H(s) are generalized non-negative functions such that 



dsG(s) < oo 



dsH(s) < oo . 



(3.3) 



Thus, we do not exclude such functions as G = 5(s — s ), < Sq < 1, etc. We shall 
see below that, for isotropic solutions (|3.1j) . each of the three equations (|2.8|h (|2.10|) . 
(I2.12jl is a particular case of Eq. (|3.2jl . 

Let us first consider Eq. (|2.8j) with f(k,t) = ip(x,t) in the notation (|3.1j) . In that 
case 



Ifcl 



, 1 ± (UJ • id) 



LO 



\k\ 



e S 



d-l 



d 



and the integral in Eq. (|2.8jl reads 

dw g(ujQ ■ ou)(p 



S d-i 



x- 



1 + c^o ■ u; 



f 



x- 



LOq ■ UJ 



(3.4) 
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It is easy to verify the identity 

/ duF{u-u ) = \S d ~ 2 \ T dzF{z){l- z 2 )^ , (3.5) 
Js*-* J-i 

where | ( S' rf— 2 1 denotes the "area" of the unit sphere in for d > 3 and IS* ) = 2. 

The identity ()3.5|) holds for any function F(z) provided the integral as defined in the 

right hand side of (|3.5|) exists. 

The integral (J3.4)) now reads 

l^" 2 | fjzg{z)(l-z 2 )^^x l ^) V {x 1 -^) = 



ds G(s)if(sx)ip[(l — s)x] , 

Jo 
where 

G(s) =2 d - 2 \S d ~ 2 \g(l- 23)13(1-3)}^ , d = 2,3,... . (3.6) 
Hence, in this case we obtain Eq. (|3.2jl . where 
(A) a(s) = s , 6(s) = 1 - s , H(s) = , 

(3.7) 

G(s) is given in Eq. (|3.6j) . 

Two other models (B') and (C), described by Eqs. (|2.12j) . (|2.10j) respectively, can 
be considered quite similarly. In both cases we obtain Eq. (|3.2|) . where 

(B') a(s) = s, b(s) = 1-8, c(s) = 1 - — ^s, 

(1 + m) z 

H{s) = 6G(s) , 

G(s) is given in Eq. (J3.6j) : while for the inelastic collision case 

(C) »W = ^, K.) = i-fi±^., 

if(s) = 0, G(s) = |5^ 2 |(l-s)^ . 

(3.9) 

Hence, all three models are described by Eq. (|3.2j) where a(s) < 1, 6(s) < 1, c(s) < 1 
are non-negative linear functions. One can also find in recent publications some other 
useful equations that can be reduced after Fourier or Laplace transformations to 
Eq. (13. 2 J) (see, for example, [1], [22] that correspond to the case G = 5(s — s ), 
H = 0). 



(3.8) 
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The Eq. ()3.2j) with H(s) = first appeared in its general form in the paper [8] in 
connection with models (A) and (C). The consideration of the problem of self-similar 
asymptotics for Eq. ()3.2j) in that paper made it quite clear that the most important 
properties of "physical" solutions depend very weakly on the specific functions G(s), 
a(s) and b(s). 



4. Models with multiple interactions and 
statement of the general problem 

We shall present in this section a general framework to study solutions to the type 
of problems introduced in the previous section. 

We assume without loss of generality (scaling transformations t = at, a = const.) 
that 



/ ds[G{s) + H(s)] = l (4.i; 
Jo 



in Eq. (|3.2j) . Then Eq. ()3.2|) can be considered as a particular case of the following 
equation for a function u(x, t) 

u t + u = T(u) , x > 0, t > , (4.2) 

where 

N N 



r(u) = «„rW(u) , « n = 1 , a n > , 

n=l n=l 

roc /*oo n 

T^ n \u) = dai . . . da n A n (ai, . . . , a n ) TT u(a k x), n = l,...,N. 
Jo Jo k=1 



(4.3) 



We assume that 

A n (a) = A n (ai, . . . , a n ) > , / da!... I da n A(a u . . . , a n ) = 1 , (4.4) 



where A n (a) = A n (ai, . . . , a n ) is a generalized density of a probability measure in R™ 
for any n = 1, . . . , N. We also assume that all A n (a) have a compact support, i.e., 

n 

A n (a!,...,a n ) = if J2 a l > R2 > n = l,...,iV, (4.5) 

fc=i 

for sufficiently large < i? < oo. 
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Eq. ()3.2|) is a particular case of Eq. ()4.2|) with 

N = 2, ai = f dsH(s), a 2 = [ dsG(s) 
Jo Jo 

1 f 1 

Ai(ai) = — / dsH(s)5[a 1 - c(s)} (4.6) 
"i Jo 

1 f 1 

A 2 (ai, a 2 ) = — / ds G(s)5[ai — a(s)]5[a 2 — b(s)\ . 

"2 Jo 

It is clear that Eq. ()4.2|) can be considered as a generalized Fourier transformed 
isotropic Maxwell model with multiple interactions provided u(0,t) = 1, the case 
N = oo in Eqs. ()4.3j) can be treated in the same way. 



The general problem we consider below can be formulated in the following way. 
We study the initial value problem 

ut + u = T(u) , u\ t =o = u (x) , x > , t > , (4.7) 

in the Banach space B = C(M + ) of continuous functions u(x) with the norm 

= sup \u(x)\ . (4.8) 

x>0 

It is usually assumed that ||mo|| < 1 and that the operator T is given by Eqs. (J4.3|) . 
On the other hand, there are just a few properties of T(u) that are essential for 
existence, uniqueness and large time asymptotics of the solution u(x,t) of the prob- 
lem (14. 7j) . Therefore, in many cases the results can be applied to more general classes 
of operators T in Eqs. ()4.7)1 and more general functional space, for example B = C(M, d ) 
(anisotropic models). That is why we study below the class ()4.3|) of operators T as 
the most important example, but simultaneously indicate which properties of T are 
relevant in each case. In particular, most of the results of Section 4-7 do not use a 
specific form (|4.3j) of T and, in fact, are valid for a more general class of operators. 

Following this way of study, we first consider the problem f)4.7|) with T given by 
Eqs. f!4.3j) and point out the most important properties of T. 

We simplify notations and omit in most of the cases below the argument x of the 
function u(x, t). The notation u(t) (instead of u(x, t)) means then the function of the 
real variable t > with values in the space B = C{R + ). 

Remark. We shall omit below the argument x e M + of functions u(x), v(x), etc., in 
all cases when this does not cause a misunderstanding. In particular, inequalities of 
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the kind \u\ < \v\ should be understood as a point-wise control in absolute value, i.e. 
< \v(x)\ for any x > 0" and so on. 

We first start by giving the following general definition for operators acting on a 
unit ball of a Banach space B denoted by 

U = {u G B : |H| < 1} (4.9) 

Definition 4.1. The operator T = T(u) is called an L-Lipschitz operator if there 
exists a linear bounded operator L : B — > B such that the inequality 

\T(u 1 )-r(u 2 )\<L(\u 1 -u 2 \) (4.10) 
holds for any pair of functions U\ j2 in U . 

Remark. Note that the L-Lipschitz condition (j4.1(Jj) holds, by definition, at any 
point x G (or x G M. d if B = C(M. d )). Thus, condition (|4.1()jl is much stronger 
than the classical Lipschitz condition 

||r(ui) -r(u 2 )|| < c\\ Ul -u 2 \\ if u li2 eu (4.ii) 

which obviously follows from (|4.10|) with the constant C = \\L\\b, the norm of the 
operator L in the space of bounded operators acting in B. In other words, the 
terminology "L-Lipschitz condition" means the point-wise Lipschitz condition with 
respect to an specific linear operator L. 

The next lemma shows that the operator T(u) defined in Eqs. (14. 3 j) . which satisfies 
T(l) = 1 (mass conservation) and maps U into itself, satisfies an L-Lipschitz condi- 
tion, where the linear operator L is the one given by the linearization of T near the 
unity. 

We assume without loss of generality that the kernels A n (ai, . . . , a n ) in Eqs. (j4.3J) 
are symmetric with respect to any permutation of the arguments (a\, . . . ,a n ), n = 
2,3,. ..,7V. 

Theorem 4.2. The operator T{u) defined in Eqs. ()4.3j) maps U into itself and satisfies 
the L-Lipschitz condition (|4.1()jl . where the linear operator L is given by 

/■oo 

Lu = / daK(a)u(ax) , (4-12) 
Jo 
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with 

N 

K { a ) = y^ria n K w (g) , 

(4.13) 

poo poo Jv 

where K n {a) — / da 2 ■ ■ ■ / da n A n (ai,a 2 , . . . ,a n ) and V^a n = l. 

Jo Jo n=1 

for symmetric kernels A n (ai, a 2 , ■ ■ ■ , a n ), n = 2, . . . . 

Proof. First, the operator T(u) in (J4.3j) - (j4.5j) maps B to itself and also satisfies 

N N 

||r(u)||<^a B ||u|| n , J>„=1. (4.14) 

n=l n=l 

Hence, 

||r(u)|| < 1 if \\u\\ < 1 , (4.15) 
and then T(U) C U, so its maps [/ into itself. 

Since T(l) = 1, we introduce the linearized operator L : B — > B such that formally 

r(l + eu) = 1 +eLn + 0(e 2 ) . (4.16) 
By using the symmetry of kernels A n (a), n — 2, 3, . . . , N, one can easily check that 

POO 

Lu = da K(a)u(ax) , 
Jo 

where 

^ poo poo 

K{a) = ^na n K n (a) , K n (a) = / da 2 . . •/ da n Ai(a, a 2 , . . . , a n ) . 

n=l J JO 

Note that K(a) > has a compact support and 

N N 



/ daK(a) = ^na„ , ^^ a n = 1 

^° n=l n=l 



because of the conditions ()4.3jl . (|4.5jl . So conditions ()4.12jl - fl4.13jl are satisfied. 

In addition, in order to proof the L-Lipschitz property 1)4.10)1 for the operator 

T given in Eqs. (|4.3|L we make use of the multi-linear structure of the integrand 

associated with the definition of T(u). Indeed, from the elementary identity 

7 7 a + c n b + d, . 
ab - cd = — — [b — a) H — (a - c) , 

we estimate 

\u\{aix)u\{a 2 x) — u 2 {a\x)u 2 {a 2 x)\ 

< \ux{axx) - u 2 (a 1 x)\ + \u 1 (a 2 x) - u 2 (a 2 x)\ , x > , a 12 > , 
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provided ||«i,2|| < 1- Then we obtain 

/■oo 

|r (2) ( Ml ) -r (2) (w 2 )| < 2 / da K 2 (a)\u 1 (ax) - u 2 (ax)\ 

Jo 

in the notation of Eqs. ()4.3|) . (|4.13j) . It remains to prove that 

POO 

|r (n) ( Ml ) - r (n) (M 2 )| < n / daiT„(a)|wi(ax) - tt 2 (ax)| (4.17) 

Jo 

for 3 < n < iV (the case n = 1 is trivial). This problem can be obviously reduced to 
an elementary inequality 

n n n 

| Y[%k - Ylvk\ < ^2 ' Xfe ~ Vk \ ' n = 3, (4.18) 

fc=i fe=i fc=i 

provided < 1, \yk\ < 1, fc = 1, . . . , n. Since this is true for n = 2, we can use the 
induction. Let 

n n 

a = x n+ i , c = y n+ i , b = x fe , d = , 

k=l k=l 

then 

n.+l n+1 

J a^/c — 2/jfc = \ab — cd\ < \a — c\ + \b — d\ < 



k=l k=l 

n 



< \x n+ i - y n+1 \ + ^2 



Xk - Uk\ , 



k=l 



and the inequality ()4.18|) is proved for any n > 3. Then we proceed exactly as in case 
n = 2 and prove the estimate ()4.18|) for arbitrary n > 3. Inequality (|4.10|) follows 
directly from the definition of operators T and L. □ 

Corollary. The Lipschitz condition (|4.11|) is fulfilled for T(u) given in Eqs. f)4.3|) with 
the constant 

N oo 

C = \\L\\ =^na n , ^2a n = l, (4.19) 

n=l n=l 

where \\L\\ is the norm of L in B. 



Proof. The proof follows directly from the inequality ()4.1(JJ) and Eqs. (|4.12|) . (J4~ 



□ 

It is also easy to prove that the L-Lipschitz condition holds in B = C(M ) for 
"gain-operators" in the Fourier transformed Boltzmann equations ([2.8)1 . ()2.9|) and 
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5. Existence and uniqueness of solutions 

The aim of this section is to state and prove, with minimal requirements, the 
existence and uniqueness results associated with the initial value problem ()4.7|) in the 
Banach space (B,\\ ■ ||), where the norm associated to B is defined in (|4.8jl In fact, 
this existence and uniqueness result is an application of the classical Picard iteration 
scheme and holds for any operator T which satisfies the usual Lipschitz condition 
(j4.11j) and transforms the unit ball U into itself. We include its proof for the sake of 
completeness. 

Lemma 5.1. (Picard Iteration scheme) If the conditions in 1)4.15)1 and 1)4.11)) are 

fulfilled then the initial value problem (J4.7)) with arbitrary u G U has a unique solution 
u(t) such that u{t) G U for any t > 0. 

Proof. Consider the integral form of Eq. (|4.2|) 

u(t)=u e- t +f dre" ( '" T) r[n(r)] (5.1) 
Jo 

and apply the standard Picard iteration scheme 

u ( n+1 \t) = n e-* + f dr e~ ( *" r) r[u (n) (r)], u (0) = u . (5.2) 
Jo 

Consider a finite interval < t < T and denote 



I m |||t — SU P ll u (^)ll • 

0<t<T 

Then 

||^ +1 )(t)||<|| M o||e-' + (l-e-')|||r( M ("))|||, 

and therefore, by induction 

||u (n) (t)|| < 1 for all n = l,2,..., and t E [0,T] , 

since ||«o|| ^ 1 an d r(tt) satisfies the inequality ()4.15j) . If, in addition, T(u) satisfies 
the Lipschitz condition (|4.11jl . then it is easy to verify that 

IK+i - "nfflr < (1 - e~ T )C\\u n - n n _i| T , n = 1,2, . . . , 



and therefore, the iteration scheme ()5.2)) converges uniformly for any < t < T 
provided 

C(l - e^ T ) < 1 =^ T < In — — if C> 1 (5.3) 

O — 1 
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(T can be taken arbitrary large if C < 1). It is easy to verify that 

u{t) = lim u {n \t) , < t < T , 



n— too 



is a solution of Eqs. ()4.7|) . f|5. 1|) . satisfying the inequality 



«(<)!!< i 



< i < T . 



(5.4) 



The Lipschitz condition (|4.11|) is also sufficient to show that this solution is unique 
in a class of functions satisfying the inequality (|5.4jl on any interval < t < s. 

This proof of existence and uniqueness for the Cauchy problem l)4.7|) - (|4.8j) is quite 
standard (see any textbook on ODEs) and therefore we omit some details. 

The important point is that in our case the length T of the initial time-interval 
does not depend on the initial conditions (see Eq. (j5.Hjl ). Therefore, we can proceed 
by considering the interval T < t < IT and so on. Thus we obtain the global in 
time solution u(t) G U, where U is the closed unit ball in B, of the Cauchy problem 



Theorem 5.2. Consider the Cauchy problem ()4.7j) with \\uq\\ < 1 and assume that 
the operator F : B — > B 

(a) maps the closed unit ball U C B to itself, and 

(b) satisfies a L-Lipschitz condition ()4.1()j) for some positive bounded linear oper- 
ator L : B -> B. 

Then 

i) there exists a unique solution u(t) of the problem ()4.7j) such that \\u(t)\\ < 1 
for any t > 0; 

ii) any two solutions u(t) and w(t) of problem ()4.7J) with initial data in the unit 
ball U satisfy the inequality 



Proof. The proof of i) follows directly from Lemma 15.11 (with the Lipschitz constant 
C = \\L\\). For the proof of ii), let u(t) and w(t) be two solutions of this problem 
such that 



□ 



u(t) - w(t)\ < exp{t(L - l)}(\u - w \) . 



(5.5) 



u(0)=u , w(0) = w , \\u \\ < 1 , ||' 
Then the function y(t) — u(t) — w(t) satisfies the equation 



w \\ < 1 . 



(5.6) 



yt + y = g(x, t) = r(u) - r(w) 



y\t=o = u -w = y . 



ON THE SELF-SIMILAR ASYMPTOTICS FOR NON-LINEAR MAXWELL MODELS 15 



Hence, 

y(t)=e-*y + [ dr e~^g(x,t) , 
Jo 

and, applying theorem 14.21 for the use of the L-Lipschitz condition (|4.1()|) , we obtain 

\y(t)\ < bole"* + f dre-^L\y{r)\ . (5.7) 
Jo 

Clearly, \y(t)\ < y*(t), where y*(t) satisfies the equation 

y.(t) = \y \e- t + [ dre-^Ly^r) . (5.8) 
Jo 

Since, by (b), the linear operator L : B — > B is positive and bounded, then Eq. (J5.8j) 
has a unique solution 

00 ±n 
n=0 

so the estimate (|5.5jl follows and the proof of the theorem is completed. □ 

Theorem 14 . 21 and the inequality (|4.14|) show that the operator V given in Eqs. ()4.3|) 
satisfies all conditions of the theorem. 

Remark. 1— The above consideration is, of course, a simple generalization of the 
proof of the usual Gronwall inequality for the scalar function y(t). The essential 
difference is, however, that y(t) is a "vector" y(x,t) with values in the Banach space 
B = C(R + ) and, consequently, the estimate (|5.5j) for the functions u(x, t) and w(x,t) 
holds at any point x G M. + . 

2— We stress that estimates ()5.7|) - (j5.8j) do not depend on specific properties of the 
operator L beyond that of being positive and bounded. The Banach space B = C(R+) 
can be also replaced, for example, by B = C{R d ) (that is the case of non-isotropic 
models). Therefore we formulated Theorem 15.21 in the most general form (without 
specifying the particular space B of continuous functions). 

We remind the reader that the initial value problem (|4.7jl appeared as a general- 
ization of the initial value problem (|3.2j) for a characteristic function tp(x,t), i.e., for 
the Fourier transform of a probability measure (see Eqs. (|3.1|) . (|2.1|) ). It is important 
therefore to show that the solution u(x, t) of the problem (J4.7)) is a characteristic 
function for any t > provided this is so for t = 0. The answer to this and similar 
questions is given in the following statement. 
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Lemma 5.3. Let U' C U C B be any closed convex subset of the unit ball U (i.e., 
u = (1 - 9)ui + 9u 2 G U' for any u lj2 G U' and 9 G [0, 1]/ If u G U' in Eq. (jUj) 
and U is replaced by U' in the condition (1) of Theorem \5. 6 A the theorem holds and 
u(t) G U' for any t > 0. 

Proof. The only important point which should be changed in the proof is the con- 
sideration of Eqs. (|5.2|) in the proof of Lemma 15.11 We need to verify that, for any 
u G U' and v(t) G U', < r < t, 

u(t) = w e"' + / dre- (i - T, !;(T) G C/' . (5.9) 

In order to see that (J5.9J) holds, we note that v(r) = T[u^(t)} in Eqs. f)5.2j) is, by 
construction, a continuous function of r G [0, t] and that 

+ [ dre^ T) = 1 . 
Jo 

Therefore we can approximate u(t) by an integral sum 

m m 

fe=i fc=i 

Then u m (t) G £/' as a convex linear combination of elements of £/'. Taking the limit 
m — > oo (f/ ; is a closed subset of U) it follows that u(t) G 17' so ()5.9|) holds. Hence, 
the corresponding sequence as in Eqs. (|52J also satisfies u^ n \t) G C/' for all n > 0. 
The rest of the proof continues as the one of Lemma 15.11 so that the result remains 
true for any closed convex subset U' C U and so does Theorem 15.21 Thus the proof 
of Lemma 15.31 is completed. □ 

Remark. It is well-known (see, for example, the textbook [17]) that the set U' C U 
of Fourier transforms of probability measures in M. d (Laplace transforms in the case 
of R + ) is convex and closed with respect to uniform convergence. On the other hand, 
it is easy to verify that the inclusion T(U') C U', where T is given in Eqs. (|4.3[) . 
holds in both cases of Fourier and Laplace transforms. Hence, all results obtained for 
Eqs. (|4.2|) . (|4.3|) can be interpreted in terms of "physical" (positive and satisfying the 
condition 1)2.1)1 ) solutions of corresponding Boltzmann-like equations with multi-linear 
structure of any order. 
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We also note that all results of this section remain valid for operators T satisfying 
conditions (|4.3J) with a more general condition such as 

N 

£}a„<l, « n >0, (5.10) 

n=l 

so that T(l) < 1 and so the mass may not be conserved. The only difference in this 
case is that the operator L satisfying conditions (j4.12j) . (|4.13|) is not a linearization of 
T(u) near the unity. Nevertheless Theorem 14.21 remains true. The inequality (|5.10j) 
is typical for Fourier (Laplace) transformed Smoluchowski-type equations where the 
total number of particles is decreasing in time (see [20,21] for related work). 

6. Large time asymptotics and self-similar solutions 

In this section we study in more detail the solutions to the initial value problem 
(I4.7|) - (|4.8J) constructed in Theorem 15.21 and, in particular, their long time behavior. 

We point out that the long time asymptotics results are just a consequence of 
some very general properties of operators F and its corresponding L, namely, that F 
maps the unit ball U of the Banach space B = C(R + ) into itself, T is an L-Lipschitz 
operator (i.e. satisfies (|4.10Jl ) and that T is invariant under dilations. 

These three properties are sufficient to study self-similar solutions and large time 
asymptotic behavior for the solution to the Cauchy problem (|4.7jl in the unit ball U 
of the Banach space C(R+). 

First of all, we show that the operator T given in Eqs. ()4.3|) has these properties. 

Main properties of the operator T: 

a) T maps the unit ball U of the Banach space B = C(R + ) into itself, that is 

\\T(u)\\ < 1 for any u G C(R+) such that \\u\\ < 1. (6.1) 

b) T is an L-Lipschitz operator (i.e. satisfies (|4.10|) ) with L from ([4.12)1 . i.e. 

POO 

\T(ui)—T(u2)\(x) < L(\ui — u 2 \)(x) = / daK(a)\ui(ax)—U2(ax)\, (6.2) 

Jo 

for K(a) > 0, for all x > and for any two functions n 12 G C(R + ) such that 

hiM\ < i. 

c) r is invariant under dilations: 

e rV T(u) = T(e TV u) , V = x^- , e rV u{x) = u(xe T ), r G R . (6.3) 
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No specific information about T beyond these three conditions will be used in 
sections 6 and 7. 

It was already shown in Theorem 15.21 that the conditions a) and b) guarantee 
existence and uniqueness of the solution u(x,t) to the initial value problem (J4.7)) - 
fl4.8j) . The property b) yields the estimate ()5.5|) that is very important for large 
time asymptotics, as we shall see below. The property c) suggests a special class of 
self-similar solutions to Eq. (14.71) . 

Next, we recall the usual meaning of the notation y = 0(x p ) (often used below): 
y = 0(x p ) if and only if there exists a positive constant C such that 

\y(x)\ < Cx p for any x > 0. (6-4) 

In order to study long time stability properties to solutions whose initial data differs 
in terms of 0(x p ), we will need some spectral properties of the linear operator L. 

Definition 6.1. Let L be the positive linear operator given in Eqs. (j4.12j) . (jHEjj) ; then 

POO 

Lx p = \{p)x p , < \{p) = / daK(a)a p < oo , p > , (6.5) 

Jo 

and the spectral function n(p) is defined by 

= M^i . (6 . 6) 

P 

An immediate consequence of properties a) and b), as stated in ()6.2|) . is that one 
can obtain a criterion for a point-wise in x estimate of the difference of two solutions 
to the initial value problem ()4.7|) yielding decay properties depending on the spectrum 
of L. 

Lemma 6.2. Let ui^{x,t) be any two classical solutions of the problem ()4.7|) with 
initial data satisfying the conditions 

\u h2 (x,0)\ < 1, Mx,0) -u 2 (x,0)\ < Cx p , x>0 (6.7) 

for some positive constant C and p. Then 

\ui(x,t) - u 2 {x,t)\ < Cx p e~ t{l - x{p)) , for all t>0 (6.8) 

Proof. The existence and uniqueness of U\^{x, t) follow from Theorem 15.21 Esti- 
mate ()5.5j) (a consequence from the L-Lipschitz condition!) yields 

lu^Xjt) — u 2 (x)\ < e~ t e Lt w(x) , with w(x) = \ui(x, 0) — u 2 (x, 0)| . (6.9) 



ON THE SELF-SIMILAR ASYMPTOTICS FOR NON-LINEAR MAXWELL MODELS 19 

The operator L from (|4.12J) is positive, and therefore monotone. Hence we obtain 

Ce x(p)t x p , 



e tL w(x) = V — L n w(x) < C e tL x p 

n 

that completes the proof. 



□ 



Corollary 1. The minimal constant C for which condition (|fi.7j) is satisfied is 



Cn 



sup 

x>0 



\ui(x, 0) - u 2 (x, 



x p 



ui(x, 0) - u 2 (x, 0)| 



and the following estimate holds 
Ui(x,t) - u 2 (x,t)\ 



< e 



-t(l-A(p)) 



•ui(x, 0) - u 2 (x, 



xp 



for any p > . 

Proof. It follows directly from Lemma ffi. 21 



(6.10) 



(6.1i; 



□ 



We note that the result similar to Lemma IB~2l was first obtained in [8] for the inelas- 
tic Boltzmann equation whose Fourier transform is given in example (C), Eq. (J2.10)) . 
Its corollary in the form similar to (|6.11|) for equation ()2.10j) was stated later in [9] and 
was interpreted there as "the contraction property of the Boltzmann operator" (note 
that the left hand side of Eq. ljOl)) can be understood as a non-expansive distance 
any between two solutions). 

However, independently of the terminology, the key reason for estimates ([6.8)1 - 
()6.11|) is the L-Lipschitz property of the operator T as defined in (|6.3j) . It is actually 
remarkable that the large time asymptotics of u(x, t), satisfying the problem ()4.7j1 
with such r, can be explicitly expressed through spectral characteristics of the linear 
operator L. 

Hence, in order to study the large time asymptotics of u(x,t) in more detail, we 
distinguish two different kinds of asymptotic behavior: 

1) convergence to stationary solutions, 

2) convergence to self-similar solutions provided the condition (c), of the main 
properties on T, is satisfied. 

The case 1) is relatively simple. Any stationary solution u(x) of the problem (|4.7jl 
satisfies the equation 



T(u) 



u . 



u E C(l 



\u\\ < 1 • 



(6.12) 
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If the stationary solution u(x) does exists (note, for example, that T(0) = and 
T(l) = 1 for T given in Eqs. ()4.3|) ) then the large time asymptotics of some classes 
of initial data Uo(x) in f)4.7|) can be studied directly on the basis of Lemma f6. 21 It is 
enough to assume that \uo(x) — u(x)\ satisfies ()6.7j) with p such that X(p) < 1. Then 
u(x, t) — > u(x) as t — > oo, for any x > 0. 

This simple consideration, however, does not answer at least two questions: 

A) What happens with u(x,t) if the inequality ()6.7|) for \uq(x) — u(x) \ is satisfied 
with such p that X(p) > 1? 

B) What happens with u(x,t) for large x (note that the estimate ()6.8|) becomes 
trivial if x — > oo). 

In order to address these questions we consider a special class of solutions of 
Eq. ()4.7|) . the so-called self-similar solutions. Indeed the property c) of T shows 
that Eq. ()4.7|) admits a class of formal solutions u s (x,t) = w(i e' 1 * 1 ) with some real 
/!*. It is convenient for our goals to use a terminology that slightly differs from the 
usual one. 

Definition 6.3. The function w(x) is called a self-similar solution associated with 
the initial value problem (14. 7|) if it satisfies the problem 

fi*Vw + w = T(w) , ||u>|| < 1 , (6.13) 

in the notation of Eqs. ()6.3|h (|4.3J) . 

Note that the convergence of solutions u(x,t) of the initial value problem (|4.7jl 
to a stationary solution u(x) can be considered as a special case of the self-similar 
asymptotics with /i* = 0. 

Under the assumption that self-similar solutions exists (the existence is proved in 
the next section), we prove the fundamental result on the convergence of solutions 
u(x,t) of the initial value problem (|4.7j) to self-similar ones (sometimes called in the 
literature self-similar stability). 

Lemma 6.4. We assume that 

i) for some /i* G R, there exists a classical ( continuously differentiable if fi* ^ 0) 
solution w{x) of Eq. (I6.13J1 such that ||u>|| < 1; 
ii) the initial data u(x, 0) = uq in the problem (|4.7jl satisfies 

u = w + 0(x p ) , ||it || < 1 , for p > such that fi(p) < /i*, (6.14) 
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where fi(p) defined in (jfi.fi j) is the spectral function associated to the operator 
L. 

Then 

\u(xe~^\t) - w{x)\ = 0{x v )e~ pt ^*-^ p)) (6.15) 

and therefore 

lim u(xe~^\t) = w(x) , x>0. (6.16) 

t— *oo 

Proof. By assumption, the function u^ix, t) = w(x e^ 1 * 1 ) satisfies Eq. (|4.7jl . Let Mi(x, i) 
be a solution of the problem (|4.7j) such that iti(a;, 0) = u (a;). Then, by Lemma IB~2l 
and by assumption ii) we obtain 

\u 1 (x,t)-w(xe fl * t )\ = Cx p e~ {1 ~ x{p))t , 

for some constant C > and all x > , i > . 

We can change in this inequality x to xe~ fJ " ,t , then 

\ Ul {xe-^\t) - w{x)\ = Cx p e-^' +l - x{p))t , 

where the tildes are omitted. Note that Ui(x,t) = u(x,t) in the formulation of the 
lemma and that p/i* + 1 — X(p) = p(/i* — fJ>{p)) in the notation ()6.6)) . 

Hence, the estimate ()6.15j) is proved. Eq. (J6.16j) follows from ()6.15|) since //* < [i{p). 
So the proof is completed. □ 

Remark. Lemma 16.41 shows how to find a domain of attraction of any self-similar 
solution provided the self-similar solution is itself known. It is remarkable that the 
domain of attraction can be expressed in terms of just the spectral function fi(p), p > 
0, defined in ()6.6j) . associated with the linear operator L for which the operator T 
satisfies the L-Lipschitz condition. 

Generally speaking, the equality ()6.16|) can be also fulfilled for some other values 
of p with /i(p) > /i* in Eq. ()6.14j) . but, at least, it always holds if (i(p) < /i*. 

We shall need some properties of the spectral function fi(p). Having in mind fur- 
ther applications, we formulate these properties in terms of the operator T given in 
Eqs. ()4.Hj) . though they depend only on K(a) in Eqs. ()6.5j) . 

Lemma 6.5. The spectral function /i(p) has the following properties: 

i) It is positive and unbounded as p — > + , with asymptotic behavior given by 

M (p) « , p^Q, (6.17) 
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where, for Y from ()4.3|) 

„oo N N 

A(0) = / da K{a) = ^ a n n > 1 , ^ a n = 1 , a„ > , (6.18) 

^° n=l n=l 

and therefore A(0) = 1 if and only if the operator V ()4.3|) is linear (N = 1); 
ii) In i/ie case o/ a multi-linear T operator, there is not more than one point 
< po < oo, where the spectral function fi(p) achieves its minimum, that is, 

dji, 
dp 1 

ClTV > 0. 



A*'(Po) — rf^(Po) — 0, with /i(po) < Mp) / or an U V > 0, provided N > 2 and 



Proof. Eqs. f!6.17|) . ()6. 18|) follow directly from the definition (|6.6|) of /x(p) and from 
Eqs. gig), (I63J) . 

The statement ii) follows, first, from the convexity of X(p) since 



\"{p)= da K(a)a p (In a) 2 >0 , (6.19) 

and from the identity 

A*'(p) = , ^(P) = pA'(p) - A(p) + 1 . (6.20) 

We note that ip{p) in ()6.20j) is a monotone increasing function of p,(ip' = p\" > 0) and 
therefore it has not more than one zero, at say, p = po > 0. Now, if TV > 2, a at > (i.e. 
T is non-linear), then p = po is also a minimum point for since from Eq. ()6.17|) 
fi(p) — > +oo as p — > and thus y!{p) < for p — »• 0. This completes the proof. □ 

Remark. From now on, we shall always assume below that the operator T from ()4.3|) 
is multi-linear. Otherwise it is easy to see that the problem ()6.13|) has no solutions 
(the condition ||w|| < 1 is important!) except for the trivial ones w = 0,1. 

The following corollaries are readily obtained from lemma 16.51 

Corollary 2. For the case of a non-linear T operator, i.e. N > 2, the spectral 
function is always monotone decreasing in the interval (0,po), and fi(p) > fi(po) 
for < p < po. This implies that there exists a unique inverse function p(/x) : 
(yu(po),+oo) — > (0,po), monotone decreasing in its domain of definition. 

Proof. It follows immediately from Lemma 16.51 part ii) and its proof. □ 



Corollary 3. There are precisely four different kinds of qualitative behavior of fi(p) 
shown on Fig.l. 
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(c) m 



Figure 1. Possible profiles of the spectral function fj,(p) 

Proof. There are two options: //(p) is a monotone decreasing function (Fig.l (a)) or 
nip) has a minimum at p = po (Fig.l (b,c,d)). In case Fig.l (a) > for all p > 
since fi(p) > 1/p. The asymptotics of A(p) ()6.5|) is clear: 

POO 

(1) A(p) >A 0O gK+ if / daiT(a) = 0; (6.21) 

p^oo J l+ 

POO 

(2) A(p) >■ oo if / daK(a)>0. (6.22) 

In the case (1) when /i(p) — > oo as p — > 0, two possible pictures (with and without 
minimum) are shown on Fig.l (b) and Fig.l (a) respectively. In case (2), from 
Eq. ()6.5|) it is clear that X(p) grows exponentially for large p, therefore fi(p) — * oo as 
p — > oo. Then the minimum always exists and we can distinguish two cases: /i(po) < 
(Fig.l (d) and /i(p ) > (Fig.l (c)) □ 
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We note that, for Maxwell models (A), (B), (C) of Boltzmann equation (Sections 2 
and 3), only cases (a) and (b) of Fig.l can be possible (actually this is the case (b)) 
since the condition (j6.21|) holds. Fig.l gives a clear graphic representation of the 
domains of attraction of self-similar solutions (Lemma 16. 4j) : it is sufficient to draw 
the line /i(p) = /i* = constant, and to consider a p such that the graph of u{p) lies 
below this line. 

Therefore, the following corollary follows directly from the properties of the spectral 
function n(p), as characterized by the behaviors in Fig.l, where we assume that 
Mpo) = for p = oo, for the case shown on Fig.l (a). 

Corollary 4. Any self-similar solution u s (x,t) = w(xe M **) with fi(po) < /i* < oo has 
a non-empty domain of attraction, where po is the unique (minimum) critical point 
of the spectral function n(p). 

Proof. We use Lemma Ifo .41 part ii) on any initial state uq = w + 0(x p ) with p > such 
that fi(po) < /i(p) < /i*. In particular, Eqs. ()6.15j) and ()6.16|) show that the domain 
of attraction of w{xe^ l * t ) contains any solution to the initial value problem ()4.7|) with 
the initial state as above. □ 

The inequalities of the kind u\ — u 2 = 0(x p ) for any p > such that fi(p) < /i*, for 
any fixed /i* > /i(po) play an important role for the self-similar stability. We can use 
specific properties of fi(p) in order to express such inequalities in a more convenient 
form. 

Lemma 6.6. For any given /i* G (/i(p ),oo) and Ui^(x) such that \\ui j2 \\ < oo, the 
following two statements are equivalent: 

i) there exists p > such that 



u 1 -u 2 = 0(x p ) , 



with ji{p) < fi* . 



(6.23) 



ii) There exists e > such that 



Ul -u 2 = 0(x v{ ^ )+e ) , with p(//„) < po , 



(6.24) 



where p(/x) is the inverse to fi(p) function, as defined in Corollary [2l 



Proof. Let property i) holds, then recall from Corollary [2] that fi(p) is monotone 
on the interval < p < po, so its inverse function p(/i) is defined uniquely. 



ON THE SELF-SIMILAR ASYMPTOTICS FOR NON-LINEAR MAXWELL MODELS 25 

It is clear, as it can be seen in from Fig.l, that the condition u(p) < /i* are satisfied 
only for some p > p(/i*), therefore inequality ()6.24j) with e = p — p(/x*) follows directly 
from ijOfl . 

Conversely, if ii) holds, first note that for any pair of uniformly bounded functions 
(note that ||wi,2|| < °o by assumption) which satisfy inequality 

\ui(x) — u 2 (x)\ < C x q , C = const. , 

for some q > 0, then the same inequality holds with any p such that < p < q and 
perhaps another constant. Therefore, if the condition (|fj.24jl is satisfied, then one can 
always find a sufficiently small < E\ < e such that taking p = p(/i*) + e\ condition 
()6.23j) is fulfilled. This completes the proof. □ 

Finally, to conclude this section, we show a general property of the initial value 
problem (|4.7|) for any non-linear V operator satisfying conditions a) and b) given in 
Eqs. (j6.1|) and ()6.2j) respectively. This property gives the control to the point-wise 
difference of any two rescaled solutions to ()4.7j) in the unit sphere of B , whose initial 
states differ by 0(x p ). It is formulated as follows. 

Lemma 6.7. Consider the problem (|4.7[) . where T satisfies the conditions (a) and 
(b) . Let Ui^ix, t) are two solutions satisfying the initial conditions Mi.2(a;, 0) = u]j 2 {x) 
such that 

||«o 2 || < 1 , «o - u l = °(0 > V > . (6.25) 

then, for any real /j,*, 

A^(x,t) = Ul {xe-»*\t) - u 2 {xe-»*\t) = 0(x p )e~ pt[ ^^ p)] (6.26) 
and therefore 

lim A„. (x,t) = , a;>0, (6.27) 

t^oo 

for any ^ > fj,(p). 

Proof. The proof is a repetition of arguments that led to Eq. (j6.15)) ). In particular, 
we obtain from Eqs. (|5.5|) . ()6.8|) the estimate 



\ Ul (x,t) -u 2 (x,t)\ = x p e t[x{p) ~ 1] ||^ 



Un ~ Ur 



and then change x to xe This leads to Eq. (J6.26|) in the notation (J6.6|) and this 
completes the proof. □ 
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Remark. There is an important point to understand here: Lemmas 16 . 21 and 16 . 41 hold 
for any operator T that satisfies just the two properties (a) and (b) stated in (j6.1j) 
and ([6.2)1 . It says that, in some sense, a distance between any two solutions with 
initial conditions satisfying Eqs. (J6.25)) tends to zero as t — > oo, i.e. non-expansive 
distance. Such terminology and corresponding distances were introduced in [18] for 
the elastic Maxwell-Bolt zmann with finite initial energy, and for specific forms of 
Maxwell-Boltzmann models in [3,22]. It should be pointed out, however, that this 
contraction property may not say much about large time asymptotics of u(x, t), unless 
the corresponding self-similar solutions are known, for which the operator T must be 
invariant under dilations (so it satisfies also property (c) as well, as stated in ([6.31) ). 
In such case one can use estimate f[6.27jl to deduce the pointwise convergence of 
characteristic functions u(x, t) (Fourier or Laplace transforms of positive distributions, 
or equivalently, probability measures) in the form ()6.15[) . (|6.16|) and the corresponding 
weak convergence in the space of probability measures [17]. We shall discuss these 
issues, in Sections 10 and 11 below, in more detail. 

Therefore one must study the problem of existence of self-similar solutions, which 
is considered in the next Section. 

7. Existence of self-similar solutions 

The goal of this Section is to develop a criterion for existence, uniqueness and 
self-similar asymptotics to the problem Eq. ()6.13)) for any operator T that satisfies 
conditions (a), (b) and (c) from Section 6, with the corresponding spectral function 
fi(p) defined in (|6.6|) . 

Theorem 17.11 below shows the criterion for existence and uniqueness of self-similar 
solutions for any operator T that satisfies just conditions (a) and (b). Then Theo- 
rem [7]21 follows, showing a general criteria to self-similar asymptotics for the problem 
()4.7j) for any operator V that satisfying conditions (a), (b) and (c) and that po > 1, 
for /i(po) = niinp>o fi{p), in order to study the initial value problem with finite energy. 

We consider Eq. (|6.13|) written in the form 



/j^xw'(x) + w(x) = g(x) , g = T(w) , /i* G M. , 



(7.1) 
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and, assuming that ||iw|| < oo, transform this equation to the integral form. It is easy 
to verify that the resulting integral equation reads 

w(x) = [ drgixr^*) (7.2) 
Jo 

We prove the following result. 

Theorem 7.1. Consider Eq. (j7.1|) with arbitrary //* G M. and the operator T satisfying 
the conditions (a) and (b) from Section 6. Assume that there exists a continuous 
function wq(x), x > 0, such that 

i) ||u> || < 1 and 

»i 

dtgo(xt^) = w (x) + 0(x p ) , g = r(w ) , (7.3) 



with some p > and p[i^ > —1 satisfying the inequality 



u(p) = - / daK(a)a p - 1 < a* . (7.4) 
P l Jo J 

Then, there exists a classical solution w(x) of Eq. (|7.1jl . 

The solution is unique in the class of continuous functions satisfying conditions 

||w|| < 1 , w(x) = w (x) +0(x pi ) , (7.5) 
with any p\ such that \i{p\) < /i*. 

Proof. The existence is proven by the following iteration procedure. We choose an 
initial approximation wq G U such that ||it»o|| < 1 an d consider the iteration scheme 

w n+ i(x) = / drg n (xr^) , p n = r(w n ) , n = 0, 1, . . . . (7.6) 



Then, property (a) of T, \\w n \\ < 1 for all n > 1 and 

|w„+i(x) - w„(a;)| < / \g n (xT^*) - g n -i(xr^)\ , n>l 



By using the inequality ()6.2|) (i.e. property (b) of T), we control the right hand side 
of the above inequality by, recalling the definition of the linear operator L from (j4.12|) 

\w n+1 {x) - w n (x)\ < / dr L(\w n - w n ^ 1 \)(xr^) 



/"l p oo 

dr da K(a)\w n — Wn-iKaxr 11 *) . (7.7) 
Jo Jo 
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Next, by assumption ii), initially 

Wi = w + 0(x p ) , with pfi* > — 1 , p > , (7.8) 

or, equivalently, 

\wi(x) — wo(x) \ < Cx p , x > , C = const. , 

then, recalling the definition for X(p) given in (|6.5jl . we can control the right hand 
side of by 

x p / daK(a)a p [ r^dr = x p A ^ . 
Jo Jo 1 + P/ 1 * 

Therefore, we estimate the left hand side of (|7.7|) by 

\w n+ i(x) - w n (x)\ < C^ n (p,n*)x p , nr(p,ti*) = z ■ (7.9) 

1 

Then, from condition ii) j(p, /i*) > since p//* > —1. Also, from condition ii), 
(Q, /i(p) = < implies < j(p,fi*) < 1. 

Therefore, there exists a point-wise limit 

w(x) = lim w n (x) (7-10) 

satisfying the inequality 

OO q 

){x) - W (x)\ < \w n+ i - W n (x)\ < r X P . (7.11) 

^—f 1-7(50, W* 



\w( 

n=0 

Estimate (|7.9|) with 7 < 1 shows that the convergence w n (x) — > w(x) is uniform 
on any interval < x < R, for any R > 0. Therefore w(x) is a continuous function, 
moreover ||w|| < 1 since \\w n \\ < 1 for all n = 0, 1, . . .. 

The next step is to prove that the limit function w(x) from ()7.1()|) satisfies Eqs. (|7.1|) . 
or equivalently, ()7.2j) . We note that 



\g n +i{x) - g n {x)\ < I daK(a)\w n+1 (ax) - w n (ax)\ < 

Jo I'-l^J 

<C\(ph n (p^*)x p . 

Therefore g n {x) — * g{x), where g(x) G C(R + ) and \\g\\ < 1, since ||g n || < 1 for 
all n. In addition, from the continuity of the operator T(u) for ||n|| < 1 follows that 
g = r(u>), and the transition to the limit in the right hand side of Eq. ()7.2|) is justified 
since \\g n \\ < 1- Hence, w(x) satisfies Eq. (|7.2|) . 

When /i* 7^ one also needs to check that Eq. ()7.1|) is satisfied. We note that, for 
any continuous and bounded w (x), all functions w n (x), n > 1, are differentiable for 
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x > and their derivatives w' n (x) satisfy the equations (see Eqs. (|7.1jl . (|7.fij) ) 

/u* x u^(a;) = iw n (x) + j^-ifa) , n = 1, 2, . . . . 

Hence, 

/i* x (w' n+1 - w' n ) = (w n - w n+ i) + (g n - g n -i) , 
and, by using inequalities (|7.9jl . (|7.12|) . we obtain 

|u4 +1 -«4| < C , 7 p ~ 1 (p,/i,) (n/(p,(J.*) + X(p))x p , n> 1 . 

Therefore the sequence of derivatives {u^x), n = 1, . . .} converges uniformly on 
any interval e < x < R. Hence, the limit function w(x) from (J7.10|) is different iable 
for x > 0, 

w'{x) = 0(x p ~ l ) , p>0, (7.13) 

and the equality (j7.1j) is also satisfied for /i* ^ 0. 

Finally we note that the condition of convergence < j(p, /i*) < 1, which is 
equivalent to the condition fi(p) < /x*, (see Eq. (J6.6|0 that has already appeared in 
Lemma 16.41 

It remains to prove the statement concerning the uniqueness of the limit function 
w(x) from (|7.1Uj) . If there are two solutions w 1 ' 2 satisfying (|7.1|) . then the integral 
equation (|7.2|) yields 

\w 1 (x)-w 2 (x)\< f drlg^xr^) - g 2 (xr' M *)\ , g 1 ' 2 = T(w 1 - 2 ) . (7.14) 
Jo 

Since ||w 1,2 || < 1, we obtain 

{w 1 (x) — w 2 (x) I < C x q , g = min(p, pi) , 

where obviously /i(q) < /i*. Then we again apply the inequality (|fj.2|) to the integral 
in Eq. (|7.2|l and get the new estimate 

K(x) - w 2 (x)\ < dx q , Ci = l{q,H*)C < C . 

By repeating the same considerations as in the existence argument, it follows that 

{w^x) — w 2 (x)\ < Cj n (q, /i*)x 9 , with j(q,fi*)<l, 

for any integer n > 0. Therefore w 1 (x) = u> 2 (x) and the proof is complete. □ 

Now we can combine Lemma f5. II with Lemma T6 .41 and prove the general statement 
related to the self-similar asymptotics for the problem (|4.7I) . 
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Theorem 7.2. Let u(x,t) be a solution of the problem ()4.7)1 with \\uq\\ < 1 and T 
satisfying the conditions (a), (b), (c) from Section 6. Let fi(p) denote the spectral 
function (17. 4|) having its minimum (infimum) at p = p (see Fig.l), the case po = oo 
is also included. We assume that there exists p £ (0,po) o,nd < e < po —p such that 

[ drgo(xr^) = u (x) + 0{x p+£ ) , g = T(u ), s>0. (7.15) 
J o 

Then 

i) there exists a unique solution w(x) of the equation (|7.1|) with /x* = fi(p) such 
that 

HI < 1 , w(x) = u (x) + 0(x p+£ ) , (7.16) 

ii) 

lim u(xe-^ p)t ,t) = w(x) , x > , (7.17) 
where the convergence is uniform on any bounded interval in M + and 

u(xe-^ p)t ,t) - w(x) = 0{x p+£ e-P {p ' e)t ) , (7.18) 

with (3{p,e) = {p + e)(fi(p) - fi(p + e) ) > 0. 

Proof. If the condition (|7.15j) is satisfied, we can take wq = Uq{x) in the assumption 
ii) of Theorem 17.11 Indeed, the function fi(p) is monotone decreasing for p £ (0,po) 
and therefore invertible. We denote the inverse function by p(/j,) and apply Lemma IB~6l 
to the condition ii) of Theorem 17.11 Thus we obtain for any /i* £ (/i(p ),oo) and 
w = u 

[ dr go(xT»*) =u + 0{x p ^" )+£ ) , e > . (7.19) 
Jo 

We note that the condition ii) cannot be fulfilled for any /i* < fi(po) since the 
set fi(p) < /i* is empty in such case. Rewriting Eq. (|7.19j) in the equivalent form 
with /i* = fi(p) we obtain the condition ()7.16|) . Then the statement i) follows from 
Theorem 17.11 Therefore we can apply Lemma 16.41 and obtain the limiting equality 
()7.17|) and the estimate ()7.18|) . This completes the proof. □ 

Thus we obtain a general criterion ()7.16|) of the self-similar asymptotics of u(x,t) 
for a given initial condition uq{x). The criterion can be applied to the problem 
()4.7j) with any operator V satisfying conditions (a), (b), (c) from Section 6. The 
specific class f!4.3|) of operators Y is studied in Section 8. We shall see below that the 
condition ()7.16|) can be essentially simplified for such operators. 
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Remark. One may expect a probabilistic connection between Theorems 17.11 and 17.21 
on the relation of rates of convergence between rates of relaxation to these selfsimilar 
states and the convergence of the corresponding Wild convolution sums formulation 
for a non-conservative Maxwell model type. Such consideration may give place to 
results of a central limit theorem type for non classical statistical equilibrium station- 
ary states. In the classical case of Maxwell Molecules models for energy conservation 
(that is for /i(l) = 0) with bounded initial energy, these connections has been fully es- 
tablished recently by [11,12], where in this case the asymptotic states are Maxwellian 
equilibrium distributions. 



The goal of this Section is to apply the general theory (in particular, Theorem I7.2jl 
to the particular case of the multi-linear operators T considered in Section 4, where 
their corresponding spectral function n(p) satisfies (J7.4)) . ()4.1Hj) and its behavior 
corresponds to Fig.l. We also show that po > 1 is a necessary condition for self- 
similar asymptotics, for p being the unique minimum of the spectral function, that 
is fi(p ) = min p>0 /i(p). 

In addition, Theorem 18 .HI establishes sufficient conditions for which self-similar 
solutions of problem (j7.1|) will lead to well defined self-similar solutions (distribution 
functions) of the original problem after taking the inverse Fourier transform. 

We consider the integral equation ()7.2|) written as 



First we establish two properties of w(x) that are independent of the specific form 
flPD of r. 

Lemma 8.1. 

i- If there exist a closed subset U' dU of the unit ball U in B, as given in (|4.15l) . 
such that r M< (U') C U 1 for any //* G 1R, and for some function u> G U' the 
conditions of Theorem \7.1\ are satisfied, then w G U' , where w is constructed 
by the iterative scheme as defined in ()7.6|) . 
ii- If the conditions of Theorem \l . l\ for V are satisfied and, in addition, T(l) = 1, 



8. Properties of self-similar solutions 





(8.1) 
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then the solution w* = 1 of Eq. (|8.1jl is unique in the class of functions w(x) satisfying 
the condition 

w{x) = 1 + 0(x p ) , /x(p) < ^* • (8.2) 

Proof. The first statement follows from the iteration scheme (|7.fi|) with u>o G t/"'. 
Then, by assumption i, wq G U' for all integer n > 1, and u> n — > w G £/'. The second 
statement follows from the obvious fact that u>* = 1 satisfies Eq. (j8.1|) with any 
/i* G M provided T(l) = 1 and from the uniqueness of w(x) stated in Theorem 17.11 
This completes the proof. □ 

The statement ii can be interpreted as a necessary condition for existence of non- 
trivial (w 7^ const.) solutions of Eq. (j8.1j) : 

if there exists a non-trivial solution w(x) of Eq. (18.1)1 for any \x*, where T(l) = \, 
such that 

|| w || = 1 , w = l + 0(x p ), p>0, then pi* < /x(p) . (8.3) 



We recall that /i(p) satisfies the inequality fi(p) > fi(po) = min p>0 /i(p) (see Fig.l). 

Moreover, if p > p (provided p < oo) in Eqs. ()8.3|) . then all solutions of Eq. (|8.1j) 
are trivial, for any /x* > fi(po), since (i(p) is increasing in (p , oo) . 

On the other hand, possible solutions with fi < fi(po) (even if they exist) are 
irrelevant for the problem (|4.7jl since they have an empty domain of attraction 
(Lemma I6.4jl . 

Therefore we always assume below that \i > n(po) and, consequently, p G (0,po) i n 
Eq. dHSD- 

Let us consider now the specific class fl4.3)) - (j4.4j) of operators T, with functions 
u(x) satisfying the condition u(0) = 1. That is, u(0,t) = 1 for the solution u(x,t) of 
the problem (|4.7J) . 

Since the operators ()4.3|) are invariant under dilation transformations ()6.3|) (prop- 
erty (c), Section 6), the problem (|4.7jl with the initial condition uq(x) satisfying 

u(0) = 1 , ||« || = 1; n (x) = l-(3x p + --- , (8.4) 

can be always reduced to the case (5 = 1 by the transformation x' = x/3 1//p . 

Moreover, the whole class of operators (|4.3|h with different kernels A n (ai, . . . , a n ), 
n = 1,2, . . ., is invariant under transformations x = x p , p > 0. The result of such 
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transformation acting on T is another operator T of the same class ()4.3j) with kernels 
A n (a 1} . . .,a n ). 

Therefore, we fix the initial condition ()8.4|) with (3 = 1 and transform the function 
(|8.4jl and the Eq. (|4.7|) to new variables x = x p . Then, we omit the tildes and reduce 
the problem (j4.7|L with initial condition (|8.4jl to the case (3 = 1, p = 1. We study this 
case, which correspond to finite energy, in detail and formulate afterward the results 
in terms of initial variables. 

Next, our goal now is to apply the general theory (in particular, Theorem 17.21 and 
criterion (|7.15Jl ) to the particular case where the initial data u (x) satisfies, for small 

x 

||u || = l, uo(x) = 1 -x + 0(x 1+£ ) , x^O, (8.5) 

with some e > 0. 

We also assume that the spectral function fi(p) given by (|7.4jl . (|4.13jl . which corre- 
sponds to one of the four cases shown on Fig.l with a unique minimum achieved at 
Po > 1- 

We shall prove that any solution constructed by the iteration scheme of theorem 
(jZHJ) with = /x(l), satisfying also the asymptotics from theorem (|7.2J) is controlled 
from below by e~ x and that lim x ^ 00 w(x) = 0, as well as other properties with a 
significant meaning. 

To this end, let us take a typical function uq = e~ x satisfying (|8.5jl and apply the 
criterion ([7. 15)) . from Theorem 17.21 or. equivalently, look for such p > for which 
(I7.15|) is satisfied. That is, find possible values of p > such that 

T Kp) (e~*)-e- x = 0(x^ £ ), (8.6) 

in the notation of Eq. ()8.1jl . 

It is important to observe that now the spectral function fi(p) is closely connected 
with the operator T (see Eqs. (|7.4jl and (|4.13|) ). since this was not assumed in the 
general theory of Sections 4-7. This connection leads to much more specific results, 
than, for example, the general Theorems 17. 1\ I7~2l 

Then, in order to study the properties of self-similar solutions and its asymptotics 
to problem (j7.1j) for p > 1, and consequently, for 

Mp)>Mpo) >-->-l, (8.7) 
Po 
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we first investigate the structure of Y^{e~ x ) for any //>—!. Its explicit formula reads 



N „ n 

r ^( e_a: ) = y2 a n / • • -dfinAnia^ . . . ,a n )Iy, iVa* 

n=l Jk + L k=l 

where 

I„(y)= dre- yT " , ^l.pO, V«„ = l. (8.9) 

1/0 n=l 

Hence, in order to find some properties of T^(y), for any \x > — 1, we prove the 
following lemma. 

Lemma 8.2. If fi > —1, y > 0, then < I^y) < 1, and 



where < r^{y) < B(y,/i) with 

(y 2 (2^+ I)" 1 if ti>-\, 

B(y,fi)= hy 2 {-hxy + y) i//i = -|; (8.11) 

V-^Py^ <f a* e (-1,-1) • 

Proof. We consider the integral I^(y), integrate twice by parts and obtain Eq. ()8.1U|) 
with 

r M (j/)=j/ 2 / e-^V^r. 



If /i > —1/2 then the estimate (j8.11|) is obvious. Otherwise we transform r y (y) into 
the form 

7>(y) = rT2/ w / e"V"H rf r ; ^ < _ 

\t L \ Jy " 

so that the estimate ()8.11|) with ji e ( — 1, — |) is also clear. 
Finally, in the case n = — ~ we obtain 

r_ 1/2 (y) = 2y 2 E(y) , with = / dr . 

Hence, rewriting 

Jl T Jy T Jy T 

C v dr 
= -\ny+ / (l- e -)- + C, 
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where C is the well-known integral 



C= / dr+ / ( e - T -l) — = / e- T lnrdr = -7, 

Ji r Jo T Jo 

with 7 = r'(l) ~ 0.577 is the Euler constant. Therefore, since C < 0, and the 

function (1 — e~ T )r~ 1 takes the value 1 at the origin and is positive decreasing in 

(0, y) for any y > 0, then the term E(y) is estimated by 

E(y) <-\ny + y. 

Hence, we obtain the estimate (j8.11|) for /1 = —1/2 and the proof is completed. □ 

Hence, we can characterize now the possible values of p > for which criterion 
(17. 15(1 holds, so Theorem 17.21 yields the self-similar asymptotics. We state and prove 
this characterization in the following corollary. 

Corollary 5. Whenever po > 1 for p,(po) = min p>0 fi(p), the condition (|8.fj|) is fulfilled 
if and only if p < 1 and, therefore, fi(p) > . 

Proof. From the previous Lemma we obtain 

I^y) = l--y— + 0(y l+£ ) , s>0, y^O, 
provided \i > — 1. Therefore 

T^e- x )-e- x = 9(fi)x + 0(x 1+£ ) , 

where 

j AT „ n 

#(//) = 1 - — — y^a ra / dai . . . da n A„(ai, . . . , a n ) V] a fe . 

1 + ^ „.=i M fc=1 

We recall that kernels yl n (ai, . . . , a n ), n = 1, . . . , N, are assumed to be symmetric 

functions of their arguments. Then 

1 f°° 
9U) = 1 A(l) , \(p) = / daK(a)aP , (8.12) 

1 + A* io 
where if (a) is given in Eqs. (|4.13|) . Recalling the definition of fi(p) in (|6.6|) . we obtain 

WO - e"* = + 0(x 1+e ) , (8.13) 

where 
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It follows from Eqs. JHHSJ), @EH1) that the condition (jSED is fulfilled if and only if 
p < 1. On the other hand, it was assumed above that po > 1, therefore > 
for such values of p. Thus, Corollary El is proved. 

□ 

In addition, by Lemma l8~Tl [ii], if /x* > //(l) then the iteration scheme 

w n+ i = T Mi (w n ) , w = e" x , n = 0, 1, . . . , (8.15) 

converges to the trivial solution it; = 1. Hence, the only nontrivial case corresponds 
to = /i(l) in Eqs. (l8~T5l . 

Hence, according to Theorem I7.1[ u> n (x) — > w(x), where u>(x) is continuously 
differentiable function on [0, oo). 

On the other hand, 

w(x) = w (x) + 0(x 1+e ) (8.16) 

and therefore w'(0) = w' (0) = —1. Since this condition is not fulfilled for w(x) = 1, 
we indeed obtain a non-trivial solution of Eq. ()8.1|) with fi* = /^(l). 

Even though we would obtain the same result starting the iterations (j8.15j) from 
any initial function satisfying Eqs. (|8.5jl . the specific function wq = e~ x has, however, 
some advantages since it gives some additional information about the properties of 
w(x). 

From now on we assume that /i* = /i(l) in Eqs. ()8.15j) and denote 

w(x) = lim w n (x) , x > . (8-17) 

n— >oo 

Then, by Theorem 17.11 w G Ci(lR + ) and satisfies the equation 

fi*xw'(x) +w(x) = T(x) , /i* = /i(l) . (8.18) 

The differentiability of w(x) was proved in Theorem 17. II only for /i ^ 0, but the proof 
can be easily extended to the case /i = since w = e~ x in Eqs. (I8.15j) has a bounded 
and continuous derivative. 

From ()8.15|) and ()8.16|) it is clear that the limit function w satisfies 

< w(x) < 1 , w(0) = 1 , w'(0) = -1 ; (8.19) 

and, by considering a sequence of derivatives in Eqs. (|8.15|) . it is easy to see that 

w'(x)<0, \w'{x)\<l. (8.20) 
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Then, estimates from Theorem 17.11 and Lemma [8.21 yield that 

w ( x ) = e ~ x + Oix^) , (8.21) 

where 

{2 for u.* > — \ , 

2 - e with any e > for = — ~ , (8.22) 
\t\ for - 1< //„<-§ . 

Hence, we collect all essential properties of w(x) in the following statement. 

Theorem 8.3. The limiting function w(x) constructed in ()8.15|) satisfies Eq. (J8-H) 
with fi = /i(l) and Eqs. ()8.18|) . where V is given in Eqs. ([4.3)1 . fi(p) is defined in 
Eqs. fl4~m TTie conditions (jOTty . (jOty . (^TTT]) are fulfilled for w(x) . More- 

over 

1 > > e~ x , lim = , (8.23) 

x— >oo 

and t/iere exists a generalized non-negative function R{t), t > 0, such that 



/•oo poo roc 

w(x)= dTR(r)e~ TX , / dr R(t) = dr R(t)t = 1 . (8.24) 
Jo Jo Jo 



Proof. It remains to prove ()8.23|) and (|8.24[) . In fact condition ()8.23|) means that e~ x 
is a barrier function to the solutions of problem ()8.1|) with [i = /i(l). 

First we note that Eq. (|8.1j) is obtained as the integral form of Eq. (17. lj) . Then, 
the identity 

fixv'(x) + -u(x) = r(u) + A(x) , 

where 

A(z) =fj,xv'(x) +v(x) -T(v) , (8.25) 
is fulfilled for any function u(a;), and the integral form of this identity reads 

v(x) = V^v) + / A(xt^ 1 ) dr . 
Jo 

Hence, if A(x) < then v < T^(v ) and vice-versa. 

We intend to prove that A(x) < for v = e~ x . If so, then w n+ i(x) > w n (x) at any 
x > in the sequence ()8.15|) generated by the corresponding iteration scheme with 
w = e~ x , and obviously w(x) > e~ x . 

Indeed, by substituting v = e~ x in Eqs. ([8.25J1 we obtain, for fi = 



= y^°nA n (x) 



n=l 
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where using (18.9jl . 

A n (x) = / da 1 ,...,da n A n (a 1 ,...,a n )Plx,'S~]ak), with 

V fe=1 / 

P(x, s) = e~ x [l - (s - l)z] - e~ sx < . 

We note that P(x, s) < for any real s and x, since e 5 ' < 1 + y for any real y. Then 
A„(x) < 0, and so also A(x) < 0. Hence, the inequality in ()8.23j) is proved. 
In order to prove the limiting identity ()8.23|) we denote 

Wno = lim w(x) . 



Such limit exists since w{x) is a monotone function. From theorem 18.31 the nice 
properties of w(x) allow to take the limit in both sides of Eq. (|8.1j) . Then 

N N 

Woo = a nw2o , ^a n = l , a n > , 

n=l n=l 

and therefore we obtain 

oo 

^a n Woo(l -w^" 1 ) = . 

n=2 

This equation has just two non- negative roots: = and = 1. The root 
= 1 is possible only if w(x) = 1 for all real x. Since by (|8.16j) this is not the case, 
then Woo = 0. 

It remains to prove the integral representation ()8.24)1 . In order to do this we denote 
by U' the set of Laplace transforms of probability measures in M + , i.e., u £ U' if there 
exists a generalized function F{r) > such that 

u{x)= drF(T)e- XT , / drF(r) = 1 . 
Jo Jo 

Then e~ x G C/' (with F = 8(t—1)) and it is easy to check that T^^U') C U' for any real 
fi. On the other hand, the set U' is closed with respect to uniform convergence in R + 
(see, for example, [17]). Thus, according to Lemma IHTT1 [i] . w G U' . On the other hand, 
it is already known from (|8.19J) that w'(0) = — 1. Hence, the corresponding function 
R{t) has a unit first moment [17]. This completes the proof of Theorem 18.31 □ 



The integral representation (J8.24)) is important for the properties of the corre- 
sponding distribution functions satisfying Boltzmann-type equations. Now it is easy 
to return to initial variables with uq given in Eq. (|8.4|) and to describe the complete 
picture of the self-similar relaxation for the problem (|4.7jl . 
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9. Main results for Fourier transformed Maxwell models 
with multiple interactions 

We consider the Cauchy problem (|4.7|) with a fixed operator T ()4.3|) and study the 
time evolution of uq(x) satisfying the conditions 

Huoll = 1 ; Uo = 1 - x p + 0(x p+£ ) , x^O , (9.1) 

with some positive p and e. Then there exists a unique classical solution u(x,t) of 
the problem (|4~7|) . (I9~T| such that, for all t > 0, 

\\u(-,t)\\ = 1 ; u(z,i) = l + 0(x p ) , x -> . (9.2) 

We explain below the simplest way to analyze this solution, in particular in the case 
of self-similar asymptotics. 

Step 1. Consider the linearized operator L given in Eqs. (|4.12|) - (|4.13|l and construct 
the spectral function given in Eq. (|7.4|) . The resulting will be of one of four 
kinds described qualitatively on Fig.l. 

Step 2. Find the value po > where the minimum (infimum) of p,(p) is achieved. 
Note that po = oo just for the case described on Fig.l (a), otherwise < po < oo. 
Compare po with the value p from Eqs. (|9.1jl . 

If p < po then the problem ()4.7jl . (j9.1j) has a self-similar asymptotics (see below). 

The above consideration shows that two different cases are possible: 

(1) P > Po provided p < oo; 

(2) < p < po, that is the spectral function /i(p) in monotone decreasing for all 
< p < p . 

In case (1) a behavior of u(x, t) for large t may depend strictly on initial conditions. 
The only general conclusion that can be drawn for the initial data ()9.1|) with p > p 
is the following: 

lim M(xe _M ,t) = 1 , x>0, (9.3) 

t— »oo 

for any [i > yu(po)- This follows from Lemma \6. 71 with = u, = 1 and from the 
fact that any such function Uo(x) satisfies the condition 

n = 1 + O(x po ) 

Case (2) with < p < po in Eqs. ()9.1|) is more interesting. In this case (assume 
that p e (0,po) is fixed) there exists a unique self-similar solution 

u 8 (x,t) =ifj(xe^ p)t ) (9.4) 
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satisfying Eqs. (|8.1jl at t = 0. We again use Lemma EHI with u\ = u and w 2 = u s and 
obtain for the solution u(x,t) of the problem (|4.7J) . (|9.1jl : 

{1 if /i > /i(p) 
^(x) if At = /i(p) (9.5) 
if /i(p) > fi > fi(p + 5) , 

with sufficiently small 5 > 0. The third equality follows from the fact that 

u (x) -tp(x) = 0(x p+e ) 
and from the equality (see Eqs. (|8.23|) ) 

lim ifj(x) = . 

x—*oo 

We note that ip(x) = w(x p ), where w(x) has all properties described in Theorem 18 .31 
The equalities ()9.5j) explain the exact meaning of the approximate identity, 

u(x,t) « ^(xe Mp) *) , t -> oo , xe M(p) * = const. , (9.6) 

that we call self-similar asymptotics. We collect the results in the following statement. 

Proposition 9.1. The solution u(x,t) of the problem (|4.7jl . (|9.1jl . wif/i T g'zwen m 
Eqs. (|4.3|) . satisfies either one of the following limiting identities: 

(1) E'g. ()9.3|) if p > po for the initial data ()9.1|) , 

(2) ()9.5|) provided < p < p . 

The convergence in Eqs. (|9.3jl . (|9.5jl uniform on any bounded interval < x < R, 
and 

u{xe^ p) \ t) - i){x) = 0{x p+£ )e-^ p ' £)t , /3(p, e) = (p + e) (a*(p) - /i(p + e) ), 
/or < p < po < e < po — P- 

Proof. It remains to prove the last statement. It follows in both cases from the 
estimate ()6.26|) for the remainder term in Lemma (|6.7|) . This completes the proof. □ 

It is interesting that our considerations are the same for both positive and negative 
values of /i(p). There are, however, certain differences if we want to consider the 
"pure" large time asymptotics, i.e., the limits (J9.3J) . (J9.5J) with \i = 0. Then we can 
conclude that 

(1) lim u(x, t) = 1 if p > po and /i(po) < , or < p < p and /i(p) < ; 

t— >oo 

(2) lim u(x, t) = ij){x) if < p < po and fi(p) = . 

t— »oo 
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It seems probable that u(x,t) —> for large t in all other cases, but our results, 
obtained on the basis of Lemma 16. 7\ are not sufficient to prove this. 

Remark. We mention that the self-similar asymptotics becomes more transparent 
in logarithmic variables 

y = \nx, u(x, t) = u(y, t) , ip(x, t) = $(y, t) 

Then Eq. f!9.6|) reads 

u(y,t) « 4>(y + fi(p)t) , t^oo, y + ^{p)t = const. , (9.7) 

i.e., the self-similar solutions are simply nonlinear waves (note that oo) = 1, 
ip(+co) = 0) propagating with constant velocities c p = — /i(p) to the right if c p > 
or to the left if c p < 0. If c p > then the value u(— oo,t) = 1 is transported to any 
given point y G M when t — > oo. If c p < then the profile of the wave looks more 
natural for the functions u = 1 — u, ip = 1 — ip. 

Thus, Eq. f!4.7p can be considered in some sense as the equation for nonlinear waves. 
The self-similar asymptotics (|9.7|) means a formation of the traveling wave with a 
universal profile for a broad class of initial conditions. This is a purely non-linear 
phenomenon, it is easy to see that such asymptotics cannot occur in the particular 
case (N = 1 in Eqs. (jUHJl) of the linear operator T. 



10. Distribution functions, moments and power-like tails 

We have described above the general picture of the behavior of the solutions u(x, t) 
to the problem (|4.7jl . On the other hand, Eq. (|4.7jl (in particular, its special 

case (|3.2j) ) was obtained as the Fourier transform of the kinetic equation. Therefore 
we need to study in more detail the corresponding distribution functions. 

We assume in this section that uq(x) in the problem (|4.7jl is an isotropic charac- 
teristic function of a probability measure in M. d , i.e., 

u (x) = F[f ] = [ dv f (\v\)e- ik - v , k G R d , x = \k\ 2 , (10.1) 

where /o is a generalized positive function normalized such that Uq(0) = 1 (distribu- 
tion function). Let U be a closed unit ball in the B = C(R+) as defined in (|4.9jl . 
Then, as was already mentioned at the end of Section 5, the set U' C U of isotropic 
characteristic functions is a closed convex subset of U. Moreover, T(U') C U' if Y is 
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defined in Eqs. ()4.3j) . Hence, we can apply Lemma and conclude that there exists 
a distribution function f(v,t), v e ~R d , satisfying Eq. (|2.1j) . such that 

, x = |A;| 2 , (10.2) 

for any t > 0. 

A similar conclusion can be obtain if we assume the Laplace (instead of Fourier) 
transform in Eqs. (jlU.lj) . Then there exists a distribution function f(v, t), v > 0, such 
that 

POO 

u(x,t) = £[/(-, t)]= / dvf(v,t)e~ xv , u(0,f) = l, x>0, t>0, (10.3) 

where u(x,t) is the solution of the problem ()4.7|) constructed in Theorem 15.21 and 
Lemma 15.31 

We remind the reader that the point-wise convergence u n {x) — > u(x), x > 0, 
where {u n , n — 1,2, . . .} and it are characteristic functions (or Laplace transforms) 
is sufficient for the weak convergence of the corresponding probability measures [17]. 
Hence, all results of pointwise convergence related to self-similar asymptotics can be 
easily re-formulated in corresponding terms for the related distribution functions (or, 
equivalently, probability measures). 

The approximate equation (|9.6|) in terms of distribution functions (J10.2J) reads 

f(\v\,t) ~ e-^ {p)t F p {\v\e-^ {p)t ) , t -»■ oo , \v\e-^ {p)t = const. , (10.4) 

where -F p (|i»|) is a distribution function such that for x = \k\ 2 

iP P (x) = T[F p \ , (10.5) 

with ip p given in Eq. (|9.4|) (the notation ip p is used in order to stress that ip defined 
in (|9.4|) . depends on p). The factor 1/2 in Eqs. (|10.4|) is due to the notation x — \k\ 2 . 
Similarly, for the Laplace transform (|10.3jh we obtain 

f{v,t) ~ e- pip)t %(ve-^ p)t ) , t -»• oo , ve~^ p)t = const. , (10.6) 

where 

^ p {x)=C\%] . (10.7) 

In the space of distributions, the approximate relation ~ is weak in the sense 
of distributions, i.e. the classical approximation concept of real valued expression 
obtained after integrating by test functions. 
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The positivity and some other properties of F p (|t>|) follow from the fact that if) p (x) = 
w p (x p ), where w p (x) satisfies Theorem 18.31 Hence 

/*oo /*oo /*oo 

i) p {x)= I dr R p (r)e- TX \ / dr R p (r) = / dr R p (r)r = 1 , (10.8) 
Jo Jo Jo 

where R p (t), t > 0, is a non-negative generalized function (of course, both ip p and 
R p depend on p) . 

We stress here that if uq{x) is a characteristic function given in Eq. (jlO.lj) . and 
condition (|9.1|) is fulfilled, then p < 1 in Eqs. (J9.1J) . In addition, the case p > 1 is 
impossible for the non-negative initial distribution f in Eqs. ()10.1|) since 

u'(0) = -C d [ dv f (\v\)\v\ 2 , (10.9) 
Jm d 

where Cd > is a constant factor that depends only on the space dimension for the 
problem. 

Hence, the self-similar asymptotics ()10.4j) for any initial data fo > occurs if po > 1 
(see Step 2 at the beginning of Section 9). Otherwise it occurs for p G (0,po) C (0, 1). 
Therefore, for any spectral function fi(p) (Fig.l), the approximate relation ()10.4j) 
holds for sufficiently small < p < 1. It follows from Eq. ()10.9|) that 

m 2 = / dv f (\v\)\v\ 2 < oo if p — 1 
Jm d 

and m 2 = oo if p < 1 . 

Similar conclusions can be made for the Laplace transforms ()10.3|) since in that 
case 

POO 

u'(0,t) = - dvf(v,t)v , 
Jo 

therefore the first moment of / plays the same role as the second moment in case of 
Fourier transforms. 

The positivity of -F(H) in Eqs. (|10.5|) and f) 1 . 7j) follows from the integral repre- 
sentation (jl0.8J) with p < 1. It is well-known that 

r-\e-\ k \ 2p )>0 , £~\e- x2p )>0 

for any < p < 1 (the so-called infinitely divisible distributions [17]). Thus, 
Eqs. ()10.8|) explains the connection of the self-similar solutions of generalized Maxwell 
models with infinitely divisible distributions. We can use standard formulas for the 
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inverse Fourier (Laplace) transforms and denote (d = 1,2,... is fixed) 





(10.10) 



J a—ioo 



Then the self-similar solutions F p and $ p (distribution functions) given in the right 
hand sides of Eqs. (jlO.fij) and ()10.8|) respectively, satisfy 



That is, they admit an integral representation through infinitely divisible distribu- 
tions. Note that Mi(|i>|) is the standard Maxwellian in R d . The functions N p (v) 
fllO.lOj) are studied in detail in the literature [17,19]. 

Thus, for given < p < 1, the kernel R p (t), t > 0, is the only unknown function 
that is needed to describe the distribution functions F p (\v |) and $ p (f). Therefore we 
study R p (r) in more detail. 

It was already noted in Section 8 that the general problem (|8.1j) . ()8.Hj) . with given 
< p < po, can be reduced to the case p = 1 by the transformation of variables x = x p . 
We assume therefore that such transformation is already made and consider the case 
p — 1. Then the equation for R(t) = Ri(t) can be obtained (see Eqs. ()8.24j) ) by 
applying the inverse Laplace transform to Eq. (J8.18|) . Then we obtain, with /z* = /z(l), 




(10.11) 





11(1)— tR(t) + R(t) = Z(R) = C-'iTiw)] , 



(10.12) 



where (see Eqs. f)4.3j) ) 

N N 



Z(R) = ^2 a n Z n (R) , ' s ^ j a n = 1 , a n > , 







k=l 
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Let us denote 



m.. 



/•oo 

/ (ItR(t)t s , s>0, (10.13) 
Jo 

then multiply Eq. ()10.12j) by r s and obtain after integration the following equality 

Gu(l)s + l)m s = 

N , n , s n (10.14) 

= ^ da 1 ... da n A n {a u ...,a n ) \dr\ . . . , dr n I ^ a fc r fe ) J R(tj) . 

n=l J*$ Vfe=l / i=l 

Next, we recall the notations flUEJ), (Q, tEE) 

„oo AT „ / n v 

Mp) = / <^ a K{a)a p = \^ a n / dai . . . da n A n (ai, . . . , a n ) I >^ a| J , 

Jo n=l M \ fc =i / 

and 

A(p) - 1 

Mp) = , 

P 

then Eq. (|10.12|) equation can be written in the form 

N 

(s/i(l) + X(s) - l)m, = ^ a n J n (s) , (10.15) 

n=2 

where 

= dai... da n A(at, ...,a n ) dr 1 ...,dT n g^\a 1 n, a n r n ) TT Rfa] 
Jm™ Jm™ 

(n \ s n 

fc=i ^ fc=i 



/■ n v s n 

g%\yi,...,y n )= [J^Vk] -J^vt , n = l,2,.... 



10.16) 



We note that gf 1 = for any s > and that mo = ni\ = 1 (see Eqs. (J10.8|0 . Our 
aim is to study the moments m s (|10.13j) . s > 1, on the basis of Eq. f)10.15j) . The 
approach is related to the one used in [22] for a simplified version of Eq. (|10.15j) with 
N = 2. The main results are formulated below in terms of the spectral function fi(p) 
(see Fig.l) under the assumption that po > 1. 

Proposition 10.1. 

[i] If the equation fi(s) = /i(l) has the only solution s = 1, then m s < oo for any 
s>0. 

[ii] If this equation has two solutions s = 1 and s = s* > 1, then m s < oo for 

s < s* and m s = oo for s > s*. 
[iii] m St < oo on/y i/ 7 n (s*) = m (|10.15j) for all n = 2, . . . , N. 
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Proof. The proof is based on the inequality 



N 



0<J2 ^nis) < C N {s)m 1 m s _ 1 (10.17) 

n=2 

with s > 1 and some positive constant Cjv(s). Then we obtain 

, C N (s) 

m s < — m s _i , m = m x = 1 . 

4/i(l) -/x(s)J 

In the case [i] /i(l) > (see Fig.l) for all s. The same is true in the case [ii] for 
s < s*. It is clear from Eq. ()10.13j) that m s > since -R(r) > and mi = 1. This 
means that the inequality ()10.17j) cannot be satisfied for s > s*, therefore moments of 
orders s > cannot exist. The statement [iii] follows directly from Eq. ([10.15)) . □ 

Hence, it remains to prove the inequality ()10.17j) . The proof is based on the fol- 
lowing elementary inequality. 

Lemma 10.2. In the notation of Eqs. ([10.16)) . 

0<g^\y 1 ,...,y n )<2 s -hL(y 1 ,y 2 ) + ^ k jP(y k+1 ,Y k )\ , s>l, (10.18) 

^ k=2 > 

•where the second term (sum) is absent for n = 2 and 

ip(yi, V2) = y{~ l V2 + , Y k = max(y 1 , . . . , y k ) , 

7^ = max(fc, fc s_1 ) , k = 1, . . . , n ; n = 2, 3, . . . 

Proof. If n = 2, we assume without loss of generality that < y 2 and reduce the 
problem (upper estimate) to the inequality 

Mx) = (1 + s ) s - 1 - x s - 2 s ~ 1 s(x + x 3 - 1 ) < , x = — < 1 . 

Its proof is obvious since, A(0) = 0, A'(x) < 0. The lower estimate in Eqs. (jl0.18|) is 
similarly reduced for n = 2 to the inequality 

g (0) = i-e s -(i-ey >o , e = —^—<\. 

x + y 

Then, its proof follows from the fact that g(0) = g(l) = 0, g"{0) < 0. 
We proceed by induction. It is easy to see that 



9n+l(yU • • • , Vn+l) = g ( n S) (yi, ...,y n )+ 92 ] ( Vn+1, ^ ^ k ) 

^ k=l ' 



n = 3, 
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Then the lower estimate (jl().18j) becomes obvious for any n > 2. By applying the 
upper estimate (jl0.18|) for #2(5) we obtain 

gnUiVl, • • • , Vn+l) < 9 { n S) + 2' 'sip ( y n+1 , Vk ) 

and note that ip(x, y) is an increasing function of y. Obviously 

n 

^2 - nYn ' Yn = max (?/i' • • • > Vn) > 
*!=i 

and therefore 

9n+i(yi, • • • , 2/n+i) < • • • , Vn) + 2 s s^ n tfj(y n+ i, Y n ) . 

This is precisely what is needed to get the upper estimate (jl(J.18j) by induction. This 
completes the proof of the lemma. □ 
In order to complete the proof of the inequality (jl0.17|) we substitute the estimates 
(|1U.18|) into the right hand side of Eq. (|10.15|) . Then the lower estimate (J10.17I) 
becomes obvious. The upper estimate (jl0.17|) also becomes obvious if we note that 

max(ai7i, . . . , a k r k ) < a n max(n, . . . , r k ) , a n = max(ai, . . . , a n ) , k = 1, . . . , n ; 

and 



j ^ dr x ... dr k ^ JJ R{Tj)J [max(ri, . . . , T k )} 1 



n\ drx... dr k I J R{jj) ) < nlm s , 

J0<r 1 <-<r fe <oo \ j=1 J 



Then it is trivial to show that the upper estimate (jl0.17|) holds. This completes the 
proof of Proposition 110.11 □ 

Now we can draw some conclusions concerning the moments of the distribution 
functions (jlU.l lj) . We denote 

1*00 poo 

m s (%)= / dv%(v)v s , m s (R p )= I dr R p (t)t s , 
Jo Jo 

m 2s (F p ) = ! dv F p {\v\)\v\ 2s , s > , < p < 1 , 

jR d 
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and use similar notations for N p (v) and M p (|u|) in Eqs. (jlO.llj) . Then, by formal 
integration of Eqs. fjlO. lljl . we obtain 

m s ($ p ) = m s (N p )m s/p (R p ) 
m 2s {F p ) = m 2s {M p )m s / p {R p ) , 

where M p and N p are given in Eqs. ()1().1()|) . 

First we consider the case < p < 1. It follows from general properties of infinitely 
divisible distributions that the moments m s (N p ) and m 2 <,(M p ), < p < 1, are finite 
if and only if s < p (see [17]). On the other hand, m>o(R p ) = m>i(R p ) = 1, therefore 
m s (Rp) is finite for any < s < 1. Hence, in this case m s ($ p ) and m2. s {F p ) are finite 
only for s < p. 

The remaining case p — 1 is less trivial since all moments of functions 



are finite. Therefore everything depends on moments of R\ in Eqs. (110. 12j) with p = 1. 
It remains to apply Proposition 110.11 Hence, the following statement is proved for 
the moments of the distribution functions (J10.4J) . (|10.6|) . 

Proposition 10.3. 

[i] If0<p< 1, then rri2 S (F p ) and m s ($ p ) are finite if and only if < s < p. 

[ii] If p = 1, then Proposition \10.1\ holds form. = m2 S {Fi) and for m s = m s ($i). 

Remark. Proposition 110.31 can be interpreted in other words: the distribution func- 
tions and & p (v), < p < 1, can have finite moments of all orders in the only 
case when two conditions are satisfied 

(1) p = 1, and 

(2) the equation n(s) = (see Fig.l) has the unique solution s — 1. 

In all other cases, the maximal order s of finite moments m2 S (F p ) and m s ($ p ) is 
bounded. 



Mx(M) = (4tt)^ 2 




d . 



Niiv) =S(v-l) , 



This fact means that the distribution functions F p and $ p have power-like tails. 
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11. Applications to the Boltzmann equation 

We recall the three specific Maxwell models (A), (B), (C) of the Boltzmann 
equation from section 2. Our goal in this section is to study isotropic solutions 
f(\v\,t),v E R d , of Eqs. (Q, Q, and (j23J) respectively. All three cases are con- 
sidered below from a unified point of view. First we perform the Fourier transform 
and denote 

u(x,t) =T[f(\v\,t)] = [ dvf(\v\,t)e- ik - v , x = \k\ 2 , u(0,t) = l. (11.1) 

JR d 

It was already said at the beginning of section 4 that u(x,t) satisfies (in all three 
cases) Eq. (|4.2jl . where N = 2 and all notations are given in Eqs. (j4.fi j) . ()3.2j) . ()3.6J) - 
()3.9j) . Hence, all results of our general theory are applicable to these specific models. 
In all three cases (A), (B), (C) we assume that the initial distribution function 

/(M,o) = / (M)>o, / tfo/ (M) = i, (n.2) 

and the corresponding characteristic function 

u(0,t) = u (x) = F[fo(\v\)} , x = \k\ 2 , (11.3) 

are given. Moreover, let Uq(x) be such that 

uq(x) = 1 -ax p + 0(x p+£ ), x->0, 0<p<l, (11.4) 

with some a > and e > 0. We distinguish below the initial data with finite energy 
(second moment) 

Eq= f dv\v\ 2 f (\v\)<oo (11.5) 
Jm d 

implies p = 1 in Eqs. (jll.4)) . as follows from Eq. (|1().9|) . and the in-data with infinite 
energy E = oo. 

Then p < 1 in Eqs. (|11.4|) and therefore 

mf = I dv f Q {\v\)\v\ 2q < oo (11.6) 

only for q < p < 1 (see [17,19]). The case p > 1 in Eqs. (jll.4j) is impossible for 
/o(l u l) ^ 0. Note that the coefficient a > in Eqs. (jll.4|) can always be changed to 
a = 1 by the scaling transformation x = a 1 ^' p x. Then, without loss of generality, we 
set a = 1 in (fTOj) . 

Since it is known that the operator T(u) in all three cases belongs to the class (|4.3)l . 
we can apply Theorem 18.31 and state that the self-similar solutions of Eq. ()4.2j) are 
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given by 

u s (x,t) = V{x) = w{x p ), (11.7) 

where w(x) is given in theorem 18 .31 and < p < po (the spectral function /i(p), defined 
in t|6.6j) . and its critical point p depends on the specific model.) 

According to Sections 9-10, we just need to find the spectral function nip). In 
order to do this we first define the linearized operator L = T'(l) for T(u) given in 
Eqs. ()4.3j) . f)4.fij) . One should be careful at this point since Ai(a\, g^) in Eqs. (|4.fij) is not 
symmetric and therefore Eqs. (|4.11|) cannot be used. A straight-forward computation 
leads to 

Lu{x) = [ dsG(s)(u(a(s)x) + u(b(s)x)) + f dsH{s)u{c{s)x) , (11.8) 
Jo Jo 

in the notation (|3.fij) - (j3.9jl . Then, the eigenvalue A(p) is given by 
L x p = A (p) x p which implies 

(11.9) 



A(p) = / rfsG( S ){(a( S )) p + (Ks)) !, }+ / ds#(s)(c( S )) p , 
and the spectral function ()6.6|) reads 

„ W = ^i. (n.io) 

p 

Note that the normalization ()4.1|) is assumed. 

At this point we consider the three models (A), (B), (C) separately and apply 
Eqs. dmSj) and (111.101) to each case. 

(A) Elastic Boltzmann Equation (Q in JR d , d > 2. By using Eqs. ()3.6|) . (|3.7j) . and 
(jOJ) we obtain 

A(p)= [' dsG(s)(s p + (l-s) p ), G(s) = A d g(l-2s)[s(l-s)}^, (11.11) 
Jo 

where the normalization constant is such that Eq. (|4.1|) is satisfied with H = 0. 
Then 

u(v) = - [ dsG(s)(s p + (l-s) p -l) 1 p>0. (11.12) 
Wo 

It is easy to verify that 

/i(p) > if p < 1 ; /z(p) < if p > 1 ; 



Ai(l) = 0, /x(2) = //(3) = - / dsG(s)s(l-s), (11.13) 

Jo 

p /i(p) — > 1 if p — > , /i(p) ^ if p — > 00 . 
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Hence, fj,(p) in this case is similar to the function shown on Fig.l b) with 2 < po < 3 
and such that fi(l) = 0. The self-similar asymptotics holds, therefore for all < p < 1. 
(B) Elastic Boltzmann Equation in the presence of a thermostat (12. 4j) in 
M. d , d > 2. We consider just the case of a cold thermostat with T = in Eq. (J3.2j) . 
since the general case T > can be considered after that with the help of (j2.11j) . 
Again, by using Eqs. (|3.6jh ([3.7)1 . and (|4.1|) we obtain 



«i pi 
A(p) = / dsG(s)(s p + (l-s) p 



Jo (l + m) 2/ 



G(s) = —A d g(l-2s)[s(l-s)}^, 
with the same constant as in Eq. (jll.ll)) . Then 

1 



11.14) 



M (p) = - / dsG(s)(s p + {l-s) p -6(l-(3s) p -(l + 6)) 
V Jo 
4m 

P = 77-, so , V > 0, 

(1 + m) 2 

and therefore 

^(1) = -0(3 [ dsG(s)s, 
Jo 

PHip) — > 1 if p — > /i(p) — > if p — > oo 



11.15) 



11.16) 



which again verifies that n(p) is of the same kind as in the elastic case A) and shown 
on Fig.l (b). A position of the critical point p such that //(po) = (see Fig.l (b)) 
depends on 9. It is important to distinguish two cases: 1) po > 1 and 2) po < 1. In 
case 1) any non- negative initial data (|11.2|) has the self-similar asymptotics. In case 
2) such asymptotics holds just for in-data with infinity energy satisfying Eqs. ()11.4|) 
with some p < po < 1. A simple criterion to separate the two cases follows directly 
from Fig.l (b): it is enough to check the sign of //(l). If 

//(l) = A'(l) - A(l) + 1< (11.17) 

in the notation of Eqs. (|11.14|) . then po > 1 and the self-similar asymptotics holds 
for any non-negative initial data. 

The inequality (jl 1.17)1 is equivalent to the following condition on the positive cou- 
pling constant 9 

0< 9<( ». = -^ SGW(slOg8 + (1 - S)1Og(1 - a)) . (11.18) 
ds G(s) (f3s + (l-/3s) log(l - (3s) ) 
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The right hand side of this inequality is positive and independent on the normalization 
of G(s), therefore it does not depend on 6 (see Eq. (I11.15J) . We note that a new class 
of exact self-similar solutions to Eq. (|2.4jl with finite energy was recently found in [10] 
for (3 = 1,6 = 4/3 and G(s) = const. A simple calculation of the integrals in ()11.18j) 
shows that 9* = 2 in that case, therefore the criterion ()11.17j) is fulfilled for the exact 
solutions from [10] and they are asymptotic for a wide class of initial data with finite 
energy. Similar conclusions can be made in the same way about exact positive self- 
similar solutions with infinite energy constructed in [10]. Note that the inequality 
fTT~T8|) shows the non-linear character of the self-similar asymptotics: it holds unless 
the linear term in Eq. ()2.4j) is 'too large'. 

(C) Inelastic Boltzmann Equation (glHJ) in M. d . Then Eqs. flUSH) and IjOjl lead to 



A(p)= f dsG(s)((as) p + (l-bs) p ) 
Jo 



where 

(1 + e) 2 , (l + e )(3-e) 



d-3 



G(s) = C d (l-s)- 
with such constant that Eq. (|4.1|) with H = is fulfilled. Hence 

i r 1 

/x(p) = -/ dsG(s)((a S y + (l-bs) p -l), p>0, (11.19) 
V Jo 

and therefore 

1 - e 2 f 1 
/i(l) = — J dsG(s)s, 

P Mp) — * 1 if p — >• , fj,(p) ^0 if p — > oo . 

Thus, the same considerations lead to the shape of u(p) shown in Fig.l (b). The 
inequality (jll.6|) with A(p) given in Eqs. (jll,19|) was proved in [9] (see Eqs. (4.26) of 
[9], where the notation is slightly different from ours). Hence, the inelastic Boltzmann 
equation (|2.5J) has self-similar asymptotics for any restitution coefficient < e < 1 
and any non-negative initial data. 

Hence, the spectral function /i(p) in all three cases above is such that po > 1 
provided the inequality ()11.18|) holds for the model (B). 

Therefore, according to our general theory, all 'physical' initial conditions (jll.2|) 
satisfying Eqs. (jll.4|) with any < p < 1 lead to self-similar asymptotics. Hence, the 
main properties of the solutions f(v,t) are qualitatively similar for all three models 
(A), (B) and (C), and can be described in one unified statement: Theorem 111.11 
below. 
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Before we formulate such general statement, it is worth to clarify one point related 
to a special value < p\ < 1 such that n(pi) = 0. The reader can see from Fig.l (b) 
that the unique root of this equation exists for all models (A), (B),(C) since = 
in the case (A) (energy conservation) , and /i(l) < in cases (B) and (C) (energy 
dissipation). If p = p\ in Eqs. ()11.4)) then the self-similar solution (jl 1.7)1 is simply a 
stationary solution of Eq. (J4.2)) . Thus, the time relaxation to the non-trivial (n ^ 0,1) 
stationary solution is automatically included in Theorem 111.11 as a particular case of 
self-similar asymptotics. 

Thus we consider simultaneously Eqs. (J2.2)) . ()2.4)1 . (J2.5)) . with the initial condition 
(jl 1.2)1 such that Eq. (jll.4)) is satisfied with some 0<p<l,e>0 and a = 1. We 
also assume that T = in Eq. ()2.4j) and the coupling parameter 9 > satisfies the 
condition (jll.lSj) . 

In the following Theorem 111.11 the solution f(\v \,t) is understood in each case as 
a generalized probability density in M. d and the convergence f n — > f in the sense of 
weak convergence of probability measures. 

Theorem 11.1. The following two statements hold 

[i] There exists a unique (in the class of probability measures) solution f(\v\,t) 
to each of Eqs. (12.2)1 . (12.4)) . (J2.5)) satisfying the initial condition (jl 1.2)1 . The 
solution f(\v\,t) has self-similar asymptotics in the following sense: 

For any given < p < 1 in Eqs. (111.4)) there exits a unique non-negative 
self-similar solution 

fi p \\v\,t) = e-^'i^Me-^') , (11.20) 

such that 

el^fQvle-^^t) -wi^M), (H.21) 

where fi(p) is given in Eqs. ()11.12j) . ()11.15jl . (jl 1.19)1 . respectively, for each of 
the three models. 
[ii] Except for the special case of the Maxwellian 

Fi(\v\) = M(\v\) = (47r)^ /2 e-T ! (H-22) 

for Eq. (j2.2)l with p = 1 in Eq. (111.4)1 (note that /x(l) = in this case), the 
function F p (\v\) does not have finite moments of all orders. IfO < p < 1, then 
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m q = dv F p (\v\)\v\ q < oo only for < q < p . (11.23) 

If p = 1 in the case of Eqs. (12. 4j) . (12.5)1 ; £/ien m g < oo only for < q < p*, 
where p* > 1 is the unique maximal root of the equation /z(p*) = /x(l), if^/i 
groen m it^s. (|1 1.4)1 . (jll.l9j) respectively. 

Proof. It is well known [17] that there exists a one-to-one correspondence between 
probability measures and their characteristics functions (111.1)1 . Moreover, the point- 
wise convergence of characteristic functions is equivalent to the weak convergence of 
probability measures. Therefore, the problem for f(\v\,t) can be reduced to study 
the problem ()4.7|) for u(x,t) with the corresponding operator T as in (J4.3j) . Then, 
the statement [i] of the theorem follows from the corresponding results for the initial 
value problem ()4.7)) (see, in particular, Proposition 19. 1)1 . We just need to express 
these results in terms of the distribution functions f(\v\,t) (see Eqs. ()10.4|l . (|10.6)) ). 

Concerning the statement [ii] of the theorem, the exceptional case of the Maxwellian 
(jll.23)) is clear since = 1 is the only root of the equation //(p*) = /i(l) = with 
fi(p) given in Eq. (|11.12j) . The rest of the statement [ii] follows from Proposition ll0.3l 
So the theorem is proved. □ 

In fact we also have some additional information (not included in the formulation 
of Theorem 111.1)1 about solutions f(\v\,t). We can express this information in the 
following way. 

Corollary 6. Under the same conditions of Theorem M 1 . 1\ the following two state- 
ments hold. 

[i] The rate of convergence in Eq. ()11.21)1 is characterized in terms of the corre- 
sponding characteristic functions in Proposition \9.1\ 
[ii] The function F p (\v\) admits the integral representation (jlO.llj) through infin- 
itely divisible distributions ([lO.lOjh 

Proof. It is enough to note that all results of section 9 and 10 are valid, in particular, 
for Eqs. (£2J), flU, (f2~5)l . □ 

Remark. To complete our presentation, we mention that a statement similar to 
theorem lll.il can be easily derived from general results of sections 9 and 10 in the case 
of 1 — d Maxwell models introduced in [15], [22], [1] for applications to economy models 
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(Pareto tails, etc.). The only difference in that case is that the 'kinetic' equation can 
be transformed to its canonical form (|4.2j) - (|4.3jl by the Laplace transform as discussed 
in section 10, and that the corresponding spectral function /j,(p) can have any of the 
four kind of behaviors shown in Fig.l. Therefore, the only remaining problem for any 
such Id models is to study them for their specific function fi{p), and then to apply 
Propositions 19.11 110.11 and 110.31 Thus, the general theory developed in the present 
paper is applicable to all existing mult i- dimensional isotropic Maxwell models and to 
1 d models as well. 
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